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Chapter 1 
Introduction 



It is well known that there are many formulations of classical mechanics, 
like the Lagrangian, Hamiltonian, Poisson bracket and Hamilton - Jacobi 
formulation 0I2]- Among these, the first two are the most important and 
most used formulations. The Hamilton - Jacobi formalism is an alternate 
formulation, which is not widely used and is more esoteric. Although, all 
these provide the same solutions to a given problem in an independent way, 
one needs to know about them because they provide different insights into 
the physics of the problem. Secondly, solving a particular problem might 
be easy in one formalism compared to the other. Hence, depending on the 
problem, one needs to choose a method, in which it can be solved easily. 
Similar to classical mechanics, there are various formulations of quantum 
mechanics [2]. Of these, some have been developed independently and some 
have their origins in classical mechanics. The quantum transformation theory 
along with the quantum Hamilton - Jacobi (QHJ) formalism, which is the 
subject of this thesis, belongs to the later group. The transformation theory 

1 



Chapter 1. Introduction 



2 



has been developed analogous to the classical theory, the roots of which, lie in 
the works of Dirac [H E] , Jordan jS] and Schwinger [Z| • The QH J formalism, 
in its present form, has been proposed by Leacock and Padgett [ El E] , on 
the lines of its classical counterpart, to study stationary states. 

In classical mechanics, the theory of canonical transformation and the 
Hamilton - Jacobi (HJ) theory are very well established These provided 

powerful formal techniques, to deal with dynamical systems, which led to a 
better understanding of the Hamiltonian formalism and the phase space. 
The HJ theory provides an efficient tool to study the systems describing 
periodic motion. Making use of the action - angle variables, one can find 
the frequency of the periodic motion, without first finding the full solution of 
the equation of motion explicitly. Similarly, in the QHJ formalism, one can 
find the energy eigenvalues without solving the Schrodinger equation for the 
stationary states. 

Analogous to the classical theory, the quantum momentum function (QMF) 
p is defined as dS/dx, where S is analogous to the Hamilton's characteristic 
function and will be called quantum action function. S satisfies the equation 



which in the limit h —>■ goes to the classical Hamilton - Jacobi equation pp . 
The solution of the Schrodinger equation, ip is related to S by 
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and one can see that p is the logarithmic derivative of ip, 



p = —ih—(\nip). 



(1.3) 



In the QHJ formahsm, the quantum momentum function plays a central 
role. Here, unlike the conventional wave mechanics, where one has to solve 
the Schrodinger equation, one makes use of the quantum Hamilton - Jacobi 
equation for the QMF 



which is a Riccati type equation. The method makes use of the singularity 
structure of the QMF to compute the energy eigenvalues in a simple and 
elegant fashion IHlEl- Thus, for quantum mechanical systems, which allow 
periodic motion, one can obtain the bound state energy eigenvalues by defin- 
ing the quantum action J, which is a generalization of the classical action 
variable and J satisfies the exact quantization condition, 



where C is the contour in the complex plane enclosing the n poles of the 
QMF on the real line. This quantization condition is exact and is a result of 
the oscillation theorem [TUl [TT] which states that the n*^ excited state wave 
function has n zeros. From ()1.3|1 . one can see that these n zeros correspond 
to the n poles of the QMF and are called moving poles, as they are dependent 
on the initial conditions. It can be shown from ()1.4|) that these poles have a 



p2 -ih^ = 2m{E -V{x)) 



(1.4) 




(1.5) 



Chapter 1. Introduction 4 

residue —ih. Thus, applying the Cauchy residue theorem to the integral in 
()4.H) . one obtains the above quantization condition. 

Along with these moving poles, the QMF will have other singularities 
originating from the potential which are called fixed singularities. Here, 
one makes a simplifying assumption that the point at infinity is an isolated 
singular point. This assumption is equivalent to saying that the QMF has 
finite number of moving poles in the complex domain. Evaluating the inte- 
gral in ()4.1|) . in terms of these other singular points and using the Cauchy 
residue theorem, one can obtain the energy eigenvalue expression. Due to 
the quadratic nature of the QHJ equation, all the residues at the fixed poles 
will have two values. Hence, one needs to choose the value of the residue, 
which will give the right physical behaviour. For this, one needs to define the 
QMF completely. The QMF is related to the classical momentum function 
Pc by the relation 

limp Pc, (1-6) 

which doubles up as a boundary condition for the QMF, as well as being 
the correspondence principle. Note that Pc is double valued and has to be 
defined precisely. The details of its definition are given in the next chapter. 

Bhalla et. al. - used this formalism to compute the bound state 
spectra of one dimensional exactly solvable (ES) supersymmetric potentials 
fSlE]- They have also demonstrated the applicability of the QHJ formalism 
to the potentials, which exhibit both broken and unbroken phases of SUSY 
(mini; depending on the range of the potential parameters. For a given po- 
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tential, it was shown that implementing the boundary conditions correctly, 
every value of the residue, at a given fixed pole was acceptable in a particular 
range of parameters and one can obtain the energy eigenvalues for both the 
phases of SUSY through the QHJ formalism. Extending this investigation 
to complicated potentials like Ginocchio, Natanzon etc, Jayanthi et. al. ob- 
tained the energy eigenvalues for these potentials and also obtained the wave 
functions for a few states [TH] . 

After successfully obtaining the solutions for bound state energies, it was 
natural to look at the other potential models like the quasi - exactly solvable 
(QES) models [12]. These are termed as QES, because one cannot analyti- 
cally obtain their entire energy spectrum but it is possible to obtain a few 
states analytically, when the potential parameters satisfy a specific condition 
known as the quasi - exact solvability condition. Applying QHJ formalism 
to these models it was found that unlike the case of ES models, where one 
obtains the energy expressions, one is lead to the quasi - exact solvability 
condition by the application of the quantization condition j20l 1^ . In this 
case, the property of square integrability was used to choose the right value 
of residue instead of the boundary condition ()1.6|) . due to the difficulty in its 
implementation to these models. 

In the above investigations except in ^Hl? the main emphasis was on 
obtaining the energy eigenvalues and not much attention was paid to eigen- 
functions. We have shown in , that one can easily obtain the bound state 
eigenfunctions for one dimensional exactly solvable models ^3] along with 
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the energy eigenstates through the QHJ formahsm. In addition to these, the 
eigenfunctions and eigenvalues for the QES models, discussed in [21]], were 
obtained in fIT\ . 

Having solved one dimensional ES and QES models, we have taken up 
the application of this formalism to more complicated problems like the pe- 
riodic potentials |2Sl-l2ni and the PT symmetric potentials [23 EH EH] • The 
application of QHJ formalism to these models is not a straight forward exten- 
sion of the method used for the bound state problems in one dimension. In 
particular, attention has to be paid to the choice of residues because square 
integrability condition cannot be insisted. Also the quantization condition 
holds as long as the oscillation theorem holds. For PT symmetric potentials, 
it is not clear whether a generalization of the oscillation theorem can be 
found. So we had to find an alternative, which will give the energy values, for 
which, the details are given in the later chapters. 

The plan of the thesis is as follows. In the next chapter, we present an 
introduction to the QHJ formalism. We will describe in detail, how one can 
employ the QHJ formalism to obtain the energy eigenfunctions and eigenval- 
ues, by explicitly working out the well known hydrogen atom potential. In 
chapter HI, the study of Scarf jHOl and Scarf - I ^3 ^] potentials is taken 
up. Both the potentials are interesting models because they exhibit different 
spectra depending on the range of the potential parameters. The Scarf po- 
tential exhibits band structure for one parameter range and ordinary bound 
state spectrum for another range. The other potential Scarf - I exhibits the 
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broken and unbroken phases of SUSY in different regions of its potential pa- 
rameters. For tliese two models, both the regions of the parameters and the 
corresponding energy spectra emerge naturally through the QHJ formalism. 

The next two chapters deal with periodic potentials belonging to the 
family of elliptic potentials [221-122]; namely the Lame potential, 

V{x) = j{j + l)m sio.'^(x,m) (1.7) 

with j being an integer and the associated Lame potential (ALP) j2Sl 

CIl { OC TJlj\ 

V{x) = a{a + l)m sn ^(x, m) + b{b + l)m ^ ^\ ' — f . (1.8) 

dn [x, m) 

Here, the functions sn (x, m), cn [x, m) and dn (x, m) are the elliptic functions 
with m being the elliptic modulus. The ALP is ES if a = b = j, j being an 
integer and is QES, when a ^ b with a + b,a — b being integers. Here, we 
show, how one can obtain the band edge solutions of these potentials using 
the QHJ formalism. This way of obtaining the solutions, turns out to be 
much simpler compared to the existing methods |21] and the comparison of 
the singularity structure with the ordinary ES and QES models is taken up 

[nE2]. 

Chapter IV deals with ES periodic potentials i.e. the Lame and the ES 
case of the ALP. For these potentials, we obtain the general forms of the 
band edge wave functions in terms of the potential parameters. Explicit 
expressions for the energy and wave functions are worked out for specific 
examples [HI]. In chapter V, we study the QES periodic potentials i.e. the 
ALP, when a ^ b and a, b taking non - integer values. Here, we obtain the 
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quasi exact solvability conditions along with the general forms of the wave 
functions [22] • 

Chapter VI deals with the study of PT symmetric ES and QES solvable 
potentials. In this chapter, the complex Scarf potential and the Khare - 
Mandal model are discussed The comparison of the singularity struc- 
ture of the QMF of these potentials with the ordinary ES and QES potential 
models is given. Chapter VII contains the concluding remarks along with a 
brief discussion on further applications of this work. 
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Quantum Hamilton - Jacobi Formalism 



In this chapter, we will briefly discuss the classical Hamilton - Jacobi (HJ) 
formalism, followed by a detailed description of the QHJ theory. By means 
of an example, the singularity structure of the QMF is explained and we will 
show, how one can obtain the energy eigenvalues from the exact quantization 
condition and simultaneously obtain the eigenfunctions. The example, which 
we use to elucidate the method, is the well known hydrogen atom, for which, 
we obtain both the eigenvalues and eigenfunctions. 

2.1 Classical Hamilton - Jacobi theory 

The classical HJ theory deals with continuous canonical transformations, 
which relates the old canonical coordinates {qi,Pi) at time t to the new 
canonical coordinates {Qi,Pi), which are constant and may have initial val- 
ues {qo,Po) at time to III 13 • The equations of transformations are of the 
following form 

qi = qiiQi,Pi,to) ; pi= Pi{Qi,Pi,to). (2.1) 
12 
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One can go from the old coordinates to the new coordinates with the help 
of a generating function. There are four types of generating functions and in 
our discussion, we consider a generating function F2{qi, Pi,t), with qi and 
as the independent variables, which will be useful for our present discussion. 
The new Hamiltonian K is related to the old Hamiltonian H by 

K = H + ^F2{q„Pi,t). (2.2) 
ot 

If K is equal to zero, then 

4 = § = 0, 

oPi 

such that, the new canonical coordinates are constants and ()2.2p becomes 

H + ^F2{q,,P,,t) = 0. (2.4) 

The transformation equations in terms of the generating function are 

d 

Pi = ■7^F2{qi,Pi,t) 
oqi 

Qi = ^F2{q,,P.,t). (2.5) 

Using the above equations in (j2.4|) . one obtains 

H{qi, g^Hqi, P., t), t) + g-/2{q^, P.. t) = 0, (2.6) 

which is known as the HJ equation, which constitutes a partial differential 
equation of P2(9i, Pi, t) in iji + l) variables. The solution of ()2.6|) is termed as 
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the Hamilton's principle function and depends on n constants of integration 
ai so that, we may write 

F2 = Siqi,ai,t) (2.7) 

and the new momenta Pi can be chosen to be the independent constants of 
integration ai, i. e. Pi = ai. Now the transformations in ()2.5|) can be written 
as 

d 

Pi = ^S{qi,ai,t), (2.8) 
dqi 

which is the classical momentum function and 

d 

Qi = -g^S{qi,a,,t)=(3i. (2.9) 

where (3i is a constant. At time to, these constitute 2n equations relating 
na's and nf3's to the initial qi,Pi values, using which, one can evaluate the 
constants of integration in terms of the initial coordinates. Thus, the above 
equation can be used to obtain qi,pi in terms of a,f3 and t as 

qi = qi{ai, f3i,t) ,pi = pi{ai, (5i,t). (2.10) 

Therefore, when solving the HJ equation, one not only obtains the generating 
function for the canonical transformation but also obtains the solution of the 
mechanical problem. 

When H is not an explicit function of t, then the HJ equation for S can 
be written as 

dS\ dS 
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where the first term involves only the q dependence and the second term 
involves time. Hence, one can assume a solution of the form 

S{qi,ai,t) = W{qi,ai) - ait, (2.12) 

where W is the Hamilton's characteristic function. Substituting this in fj2.1ip 
gives the HJ equation for the characteristic function as 

ff(,..g^j=a. (2,13) 

The transformation equations will be 

dW 



dqi 

Since, H and W are independent of t, from ()2.2|1 . one obtains 

K = ai (2.15) 

and the equations of motion are 

OK 

with solution Pj = and 

4 = §^ = ii = i 



with solutions 



2 ^ 1 (2.17) 



Qi = t + , 2 = 1 

dW 

Q, = = (2.18) 
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From the above equation, we see that only Qi is not a constant, but is a 
function of time. Thus, with HJ equation of the characteristic function, one 
can solve the systems in which H is independent of t. 

2.1.1 Action - angle variables 

As a variation to the above discussed HJ method, one can deal with systems 
with periodic motion. Here, for each periodic motion, one introduces the 
action variable defined in terms of the momentum p as 

J = j^pdq, (2.19) 

where the integration is carried over a complete period of the periodic motion 
in the phase space. The Hamilton's characteristic function in terms of this 
new variable is 

W = W{q,J). (2.20) 

The coordinate conjugate of J, known as the angle variable is defined as 

dW 

and the corresponding equation of motion is 

u = = V (2.22) 

with the solution 

^ = vt^^ (2.23) 

where v is the frequency. Equation ()2.2H1 can be solved for g as a function 
of ijj and J, which, in combination with ()2.23|) . will give the desired solution 
as a function of time. 



Chapter 2. Quantum Hamilton - Jacobi Formalism 17 

In order to see how these variables are advantageous, consider the change 
in uj, as q goes through a complete cycle of rotation, which is given by 




(2.24) 



using (j2.2H) . one gets 




(2.25) 



Since J is a constant, one can write 

where ()2.19|) has been put to use. From ()2.23|) . it follows that, if r is the 
period for complete cycle of q, then 

Auj = l = UT. (2.27) 

Hence, the constant z/ can be identified as the reciprocal of the time period 

y=- (2.28) 

T 

and therefore, the frequency is associated with the periodic motion of q. Thus 
for a one dimensional periodic system, the frequency of the system can be 
found by determining if as a function of J and using ()2.22|) . Hence the use 
of action - angle variables, provides a powerful technique for obtaining the 
frequency of the periodic motion, without finding a complete solution of the 
motion of the system. 
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2.2 Quantum Hamilton - Jacobi Formalism 

The QHJ formalism draws its inspiration from the classical HJ theory and 
was first initiated in this form by Leacock and Padgett, in 1983 j3l4j . As 
mentioned in the introduction, the central entity in this formalism is the 
quantum momentum function p, which is defined in three dimensions as 

p = VS (2.29) 

where S is analogous to the classical characteristic function W. S is related 
to the solution of the Schrodinger equation, 

- 7^V^?A(x, y, z) + V{x, y, z)i){x, y, z) = Eip{x, y, z) (2.30) 



i^{x,y,z) = expi — ] (2.31) 



by the relation 

lh(T.. 7.\ = PTCn I 

h 

which, when substituted in (j2.3(Jj) . gives 

(V^)^ - ihV.iVS) = 2m{E - V{x, y, z)). (2.32) 

Substituting (j2.29j) in (|2.32j) gives the QHJ equation for p as 

{pf -inV.p = 2m{E -V{x,y,z)) (2.33) 

and from (j2.29|) and (j2.31|) . one can see that p is the the logarithmic derivative 
of 'ipi^x, y, z) i. e. 

Vip{x,y,z) 

p = -ih— (2.34) 

i)[x,y,z) 
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The above discussion of the QHJ formahsm is done in three dimensions but, 
from here on, we proceed to one dimension, as we will be discussing only one 
dimensional potentials in this thesis. One can see that in one dimension, the 
QHJ equation ()2.HH|1 takes the form of the Riccati equation 

-ih^ = 2m(E -V(x)). (2.35) 

ax 

In the QHJ formalism, x is treated as a complex variable, there by extending 
the definition of p to the complex plane. To obtain the full solution of the 
Riccati equation, one needs to know p completely in the complex plane. From 
()2.33p . one observes that 

lim p — > pc, (2.36) 

where Pc = \j2m{E — V{x)) is the classical momentum, which can be used 
as a boundary condition on p. When x is a complex variable, Pc needs 
to be defined carefully, since, it is a double valued function. In the phase 
space, Pc vanishes at points where E = V{x). In the complex x plane, the 
classical region is on the real axis, in between the turning points xi and X2- 
Correspondingly, pc is given a branch cut in between the turning points and 
is defined to be the branch, which is positive just below the branch cut jHl 
H] . This completes the definition of Pc and the role of ()2.36|) is not only that 
of a boundary condition on p, but also as a correspondence principle. 

2.2.1 Singularities of the QMF 

It has already been mentioned that the singularity structure of the QMF plays 
an important role in obtaining the eigenfunctions and eigenvalues through 
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the QHJ formalism. The location of singularities of p and the residues at 
these singular points are the essential ingredients required in this formalism. 

The QMF has two types of singularities, the moving and the fixed singu- 
larities. It is well known from the oscillation theorem, that the n*^ excited 
state wave function has n zeros on the real line. These zeros correspond 
to the poles of p, as seen from ()2.34|) . The location of these poles depends 
on the energy and the initial conditions. Hence, these are called as moving 
poles. For Riccati equation, it is known that only poles can appear as moving 
singularities It is easy to see that the residue at these moving poles is 
—ih. Let xq be a point, where V{x) is analytic and p has a moving pole. 
Assuming a moving pole of order m, we proceed to calculate the residue of 
p at xq, by expanding p in Laurent expansion around xq as 



Substituting ()2.37|1 in the QHJ equation ()2.35j) and comparing the coefficients 
of different powers of x — Xq, one gets m = 1 and 



Thus, the moving poles are simple poles with residue —ih. This information 
of the moving poles is used to obtain the quantization condition satisfied by 
the quantum action 



Here, C is the contour, which encloses the moving poles of p as shown in the 





Ai = —ih. 





(2.39) 



fig. 2. L 
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Fig. 2. 1 



Applying the Cauchy residue theorem to the n poles enclosed by C, one 
obtains the quantization condition 



J{E) = — <j> pdx = nh, 



(2.40) 



which can be viewed as a generalization of the Wilson - Sommerfield quanti- 
zation. Unlike the WKB quantization and the Wilson - Sommerfield quanti- 
zation, this is an exact quantization condition and is a result of the oscillation 
theorem and is true only for the bound states OlSl- Hence one cannot ap- 
ply ()2.40|) to continuous energy states with wave functions non vanishing at 
infinity. 

The power of the classical HJ theory lies in the fact that one obtains the 
frequencies of the system directly from the classical action - angle variables. 
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Analogously, in the QHJ formalism, one can obtain the energy eigenvalues 
of the stationary states directly from the quantum action variable making 
use of the exact quantization condition ()2.40|) . To arrive at the energy eigen- 
values, one needs to know the location of other singularities of the QMF, 
which lie outside the contour C and their corresponding residues. The other 
singularities of the QMF arise from the potential V{x) and are called as fixed 
singularities. The locations of these singular points are energy independent 
and appear in all the solutions of the QHJ equation. Thus, if the potential 
is meromorphic, the QMF will also be meromorphic. The residue at any of 
these fixed poles can also be computed by using the QHJ equation ()2.H3|1 . 
Due to the quadratic nature of the QHJ equation, one obtains two values for 
each of these residues at the fixed poles. Of these two values, one needs to 
check which of the values satisfy the boundary condition ()2.36|1 . Since, Pc 
is defined to be positive just below the branch cut, we choose those values 
of residues, which in the limit h 0, correspond to positive pc below the 
branch cut. This method of choosing the residue is explained in the study of 
hydrogen atom, later in the chapter. 

In all our investigations, we have assumed that the QMF has only finite 
number of singularities in the complex plane, which turned out to be true for 
all the ES and QES models studied. This assumption is equivalent to saying 
that the point at infinity is an isolated singularity. This property is shared 
by both the ES and QES models and the only difference is the location of 
the moving poles. For ES models, the QMF has poles only on the real line, 
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whereas for QES models, the moving poles are in both real and complex 
locations. In case of periodic potentials, it will be shown that both ES and 
QES cases share the same kind of singularity structure. The knowledge of 
the singularity structure allows one to write the QMF as a function of x and 
hence compute the eigenvalues and the eigenfunctions. The integral J{E) in 
()2.40|) can be computed in terms of the residues at the singularities, outside 
the contour C. This gives the energy eigenvalues for one dimensional ES 
models [7] and for QES models, one is lead to the QES jHl E] condition by 
the same condition. 

In the next section, we clearly explain all the above discussed points using 
the example of the hydrogen atom. We also show directly from the QHJ 
formalism that without using the quantization condition, we can arrive at the 
form of the QMF and hence, the values of the eigenvalues and eigenfunctions 
by using well known results in complex variables. 



2.3 Hydrogen atom 

In this section, we obtain the bound state wave functions for the radial part 
of the Coulomb problem PU] • We consider the supersymmetric potential 

Vir) = -'- + ^l±^+ (241) 

4 

which is equal to the original coulomb potential plus the constant 
which has been added to make the ground state energy equal to zero. The 
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radial Schrodinger equation with h = 2m 



1 is given by 



(PR 2dR ( 

\ \- 




) 



R = 0. 



(2.42) 



dr'^ r dr \ r 



4(/ + l)2 



Using the transformation R{r) = u{r)/r, p = —i{u'{r)/u{r)), one obtains the 
QHJ equation as 



The wave function u{r) should vanish at r = and hence the QMF has fixed 
pole at the origin along with the n moving poles with residue —i on the 
positive real line. We assume that there are no other singular points of p in 
the complex plane, which makes the point at infinity, an isolated singularity 
of p. With this singularity structure of p, we, now proceed to find the energy 
eigenvalues. 

From this information about the singularity structure of the QMF, we 
know that p is meromorphic and has the above mentioned singular points. 
Hence, one can write p as a sum of singular and analytic parts as written 
below 



In the above equation, the summation term describes the sum of all the 
principal parts coming from all the individual Laurent expansions of p taken 
around each moving pole with —i as the residue. The second term in (I2.44j) 
represents the singular part in the Laurent expansion of p around r = 0, with 
bi as the residue and Q{r) represents the analytic part of p. From the QHJ 
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equation ()2.4Hj) . one can see that for large r, p goes to zero. Thus Q{r) is 
an entire function and bounded at infinity. Thus, from Louville's theorem, 
Q{r) is a constant, say C. Thus, fl2.44p becomes 

P = E^^ + - + C7, (2.45) 

r - rfe r 

which gives the form of p in the entire complex plane. 

To obtain the energy eigenvalues one only needs to calculate the value 
of the residue bi and the value of the constant C. On the other hand, the 
location of the moving poles must be found in order to obtain the wave 
function. 

2.3.1 Definition of pc 

As mentioned in the previous section, due to the quadratic nature of the QHJ 
equation, bi will have two values. To select the right value of residue, we make 
use of the boundary condition ()2.H(jj) . which uses the classical momentum pc- 
Hence, we first proceed to define pc which is given by 



Pc 



.,f!^i(l±i)__£!_ (2,46) 

\ r r2 4(/ + l)2 



in the complex plane and which is equivalent to 

±ia ( 2 eV /(/ + 1)^^/^ 



where a = ^ ^^^^-^^2 — E, is a positive constant. Let ri and r2 be the two 
turning points and in the complex r plane, Pc is given a branch cut between 
the two turning points. Pc is defined to be positive just below the branch cut. 
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Hence, one needs to choose that branch of pc from ()2.47j) which satisfies this 
condition. In order to be able to do that we introduce the variables, pi, p2, 
and 02 such that 

r - ri = pi exp(z6'i), r - r2 = P2 exp(z6'2) (2.48) 
where 6i and 62 are the angles as shown in the fig 2.2 below 

I(x) 



r 




r r 
1 2 



Fig. 2. 2 

and are chosen to lie in the open interval (0, 2tt). This gives a branch cut to 
Pc, for ri < r < r2. Thus, one can write Pc as 

Pc ~ ^V'(r-ri)(r-r2), (2.49) 
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where A = ±1 and using ()2.48|1 . pc becomes 



Pc = -^ipip2) ' exp \ i ^ 1 . (2.50) 

We need to fix the sign of A, which will fix the branch of Pc to be selected. 
For this, we consider a point just below the branch cut which corresponds to 
Q\ ~ 27r, Q2 ~ TT, which gives 

Pc = ^{pxP2)"^. (2.51) 

Therefore, for p^ to be positive just below the branch cut, we choose A = +1. 
Thus, we obtain the definition of pc in the entire complex plane as 

Pc = ^{Pip2) ' exp U ^ 1 , (2.52) 

with < < 277 and < ^2 < 27r. 

2.3.2 Residue at the fixed pole r = 

With the complete definition of Pc obtained, we proceed to calculate the 
residue of both Pf. and p at r = 0. Comparing the residues, one can fix the 
correct value of 61. First, we calculate the residue of Pc- At the origin ^1 ~ tt 
and Bi ~ vr, which, when substituted in fl2.52|l . give the residue at r = as 
(E3^ at r = 

'"(rirs)^/^ (2.53) 



2 

Now to calculate the residue of p, we perform a Laurent expansion of p around 
r = as 

p = — + ao + fli'^ + (2.54) 

r 
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Substituting ()2.54j) in ()2.4H|1 and comparing the coefficients of one ob- 
tains the quadratic equation 

bj + ibi + lil + l) =0, (2.55) 

which gives the two values of bi as 

bi=il , + 1). (2.56) 

Comparing the values of bi with the residue of pc in ()2.53|) . we see that the 
sign of the imaginary part must be negative. Hence, we must select 

6i = -i(/ + l), (2.57) 

which satisfies the boundary condition ()2.36|) . 

2.3.3 Calculation of energy eigenvalues and eigenfunctions 

With the selection of the right value of the residue at the fixed pole, one can 
write p as 



where J2k=i ^ri^^ = ""^p^ and P„(r) is an n*^ degree polynomial. Substi- 
tuting p from ()2.58p in ()2.43|) one obtains 

Pn{r) \Pn{r) r J r Pn{r) 
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In order to fix C and obtain the energy eigenvalues, we look at the behaviour 
of each term of ()2.59j) . for large r. Noting that for large r, 

P^ir) n{n-l) P^(r) n 
-Pn(r) Pn[r) r 

we equate the coefficients of the constant terms and coefficients of 1/r in 
fl2.59|l to zero individually and get 

C-^^^:^yE-C^^^^ (2.61) 

which in turn give 



4 4 

e e 



^- 4(^ + / + i)2 + 4(; + i)2' (2-62) 

which is the required energy eigenvalue expression and one is left with the 
second order differential equation 

Y r n + L + lj yr^ [n + l + \)rj 

(2.63) 

Defining ^q^^r = y in ()2.63|) . we get 

+ ((2/ + !) + !- y)P^(i/) + nPn{y) = (2.64) 

which is the associated Laguerre equation, and P„(y) is proportional to the 
Laguerre polynomial, denoted by L^'+^(y). The bound state wave function 
can be obtained by using ()2.34|1 . which gives 

u{r) = exp j pdx^ (2.65) 
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Substituting the value of p from ()2.58|1 in the above equation and changing 
the variable from r to y, one obtains the wave function in terms of y as 

u„(y) = y'+iexp(-|)L^i(y), (2.66) 

which matches with the known result [TT] . 

We would like to mention here that one can obtain the energy eigenvalues 
by using the exact quantization condition ()2.40|) discussed in §2.2 without 
solving the QHJ equation. For this, one needs to calculate the residue at 
infinity , which will again have two values. The right value is chosen by 
applying the boundary condition ()2.36|) and by using the Cauchy residue 
theorem, one is led to the expression for the energy eigenvalues, the details 
of which are given in ^2]- Instead of using the above mentioned boundary 
condition, one may note that p ()2.45|) is a rational function of x and hence, 
the sum of all the residues (including the residue at infinity) is zero . This 
is equivalent to what one gets from the exact quantization condition and is 
widely used in our analysis of periodic and PT symmetric potentials. 

2.4 Change of variables in QHJ Formalism 

We would like to emphasize here that, locating the singular points of the 
QMF and imposing the boundary condition ()2.36|1 are two crucial steps in 
arriving at the correct physical solutions. For many problems, it is necessary 
to perform a change of variable 

y = fix). (2.67) 
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Thus, for a general potential V{x), the QHJ equation ()2.35|1 with h = 2m = 1 
is 

p"^ -ij^-{E- V{x)) = 0, (2.68) 

which, after a change of variable ()2.fi7j) and introduction of q by q{y) = 
ip{x{y)), becomes 

q'iy) + Fiy)^ + E- Viy) = (2.69) 
ay 

Here F{y) equals the derivative expressed as a function of y and V{y) = 
V{x{y)). In order to calculate the residue at the moving poles, removing the 
complexity caused due to the presence of F{y), we bring ()2.69|) to a simple 
form, by introducing x through the following transformation equations. 

1 F'(v) 

q = F{y)<P, ;^ = 0+--^, (2.70) 

when substituted in ()2.69j) . give 

dx E-V(y) UF'(y)\ llF(y)V 

where, the residue at moving poles is one. The function x will also be called 
as the QMF and equation ()2.71|) as the QHJ equation. From ()2.71|) . one can 
again see that the residue at the moving poles of x is unity. The expression 
for the wave function, in the new variable y can be obtained as follows. From 
()2.34p . one obtains ip with % = 1 a,s 

ijj{x) = exp (^i J pdx^ (2.72) 
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in which, substituting p = —iq, followed by the change of variable ()2.(j7j] 
gives 



Using the transformation equations given in ()2.7()j] . one is led to the expres- 
sion for ipiy) ill terms of x and F{y) as 



It is not always easy to impose the boundary condition (j2.36|) . We, there- 
fore, describe a set of constraints, which will be found useful in restricting 
the choice of residues at the fixed poles. Firstly, the wave functions should 
be finite every where for the physical values of x. If one is looking for bound 
state solutions, one must have the condition that the wave function vanishes 
for large x. Since, in this thesis, we are dealing only with the one dimen- 
sional potentials, we can also use the fact that vanishing of wave function 
for large x, implies that the bound state levels are non - degenerate. So, for 
example, parity invariance will imply ■?/'(— x) = ±iIj{x) and QMF will satisfy 
p{—x) = —p{x). Also the bound state wave functions will be real and hence 
p*{x) = p{x). All possible values of the residues, consistent with these and 
other similar conditions, must be accepted. In all the cases studied here, the 
QMF turns out to be a rational function and hence, the quantization condi- 
tion becomes equivalent to the sum of all the residues of QMF, including the 
residue at infinity, being equal to zero. This leads to a relation between the 
residues and the number of moving poles n. As we will see in later chapters. 




(2.73) 




(2.74) 
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some combinations of residues get eliminated because the number n lias to 
be positive. 
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Chapter 3 



Study of potentials exhibiting different 

spectra 

3.1 Introduction 

This chapter is devoted to the study of potentials, which exhibit different 
spectra for different ranges of the potential parameters. The first potential 
investigated, here, is the Scarf - I potential, which has been well studied [H 
12] and is known to exhibit both phases of SUSY, namely the broken and the 
unbroken phases of SUSY, for different parameter ranges. In Bhalla et. 
al. have studied this potential using the QHJ formalism and obtained the 
the energy eigenvalue expressions in both the phases. In the present study, 
we obtain both the eigenvalues and eigenf unctions along with the different 
ranges. 

The second potential to be discussed is the Scarf potential|ll E], which 
exhibits band structure in one range and ordinary bound state spectrum in 
the other range. In [S], a group theoretical treatment of Scarf potential has 
been discussed with the bound state and scattering state problem, being 

36 
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extended to include the band structure. Using the QHJ formahsm, we show 
that one can arrive at the bound state solutions and the band edge solutions 
of this potential in a simple and straight forward fashion. The next section 
of the chapter discusses the Scarf - I potential, followed by a section on Scarf 
potential. 

3.2 Phases of Supersymmetry 

The expression for the supersymmetric Scarf - I potential is 

VI (x) = -A^ + {A^ + - Aah) sec^ ax - B(2A - ah) tan ax sec ax (3.1) 

and it exhibits broken and unbroken phases of SUSY for the following ranges 
of potential parameters p. In the parameter range, 

{A-B)> 0, {A + B)> 0, (3.2) 

SUSY is exact and SUSY is broken for the range 

{A-B)> 0, {A + B) < 0. (3.3) 

The supersymmetric phases of this potential are well studied and their cor- 
responding bound state solutions are available in the literature [Dl^ . In this 
chapter, we will show that the ranges of the two different phases and their 
corresponding solutions emerge naturally from the QHJ analysis [H]. 

Putting 2m = 1 and writing the QHJ equation in terms oi q = ip we get 

+ hq' + E + A^ - {A^ + 5^ - Aah) sec^ ax + B{2A - ah) tan ax sec ax = 

(3.4) 
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As discussed in the end of §2.4, we perform a change of variable 

y = sin ax, (3.5) 

which gives F(y) = ay/1 — y^. The transformation equations (j2.7Up of §2.4 
take the form, 

2(1-1/2 



q = a^l-y^(f) ; (f) = X + 7^- (3-6) 



and the QHJ equation in terms of x is obtained as, 

2 , , y'^h^ E + A^ a'^h - 2(A'^ + - Aah) 



2\2 



4(1-1/2)2 «2(1_^2) 2^2(1 

, ^?<^i^^0. (3.7) 

3.2.1 Form of QMF (x) 

From ()3.7p . one can see that x has fixed poles at ?/ = ±1. In addition, there 
will n moving poles, with residue one, on the real line corresponding to the n 
nodes of the n^^ excited state. Proceeding in the same way as in the case of 
hydrogen atom, we assume that x has only the above mentioned singularities 
in the complex plane. Hence, we write x 

, = A^ + 4^ + fi5M+c, (3.8) 

y-l y+l Pn{y) 

where bi and b'l are the residues at y = ±1 respectively and Pn{y) is an ra*'* 
degree polynomial. C represents the analytic part of x which is a constant 
due to Louville's theorem, since x is bounded at infinity as seen from ()3.7|1 . 
As in §2.3, we first define the classical momentum function, which will help 
us to fix the right form of x- 
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3.2.2 Definition of the Classical momentum function 

The classical momentum function, which is given a branch cut in between 
the two turning points Xi and X2 in the complex x plane, is given by 



= pm{E - V{x)) (3.9) 

and is defined to be the branch, which is positive just below the branch cut. 
Since all the analysis is performed in the complex y plane with x as the QMF, 
which has fixed poles at y = ±1, we use the condition 

limx^Xc (3.10) 

instead of ()2.36|) . The expression of Xc is obtained as, 

ipc yh 
Xc = / -I ■) ~ 2^' (3-11) 

aVl-r 2(1-?/^) 

by putting q = ip in the transformations ()3.(jj) and replacing p with pc, which 

represents the classical momentum function ()3.20|) in terms of y as 

_ , . VeTA^ f , , B'~Aa-E-Bi2A-Aah)y Y^' 

= ^^7r^ V ^Ta^ ) ' ^ ^ ^ 

where 1 — is positive as y lies between ±1. Let yi and y2 be the images 
of the two turning points Xi and X2 in the y plane. The mapping from 
y = sin ax in (j3.5p . is such that it maps the point just below the branch cut 
in the x plane to a point just below the branch cut in the y plane. Hence, we 
select the branch of Pc, which is positive just below the branch cut. Similar 
to the hydrogen atom case, we introduce the variables 9i,ri and 62, r2, where 

y - yi = nexpiOi y - y2 = r2exp 162 (3.13) 
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as represented in the fig. 3.1. 




Fig. 3. 1 

The values of 6i and 62 are chosen to he in the interval (0, 2tt) and pc can 
be written in terms of the above variables as 



Pc = ,, -(rir2) ' exp i 



(3.14) 



where A = ±1. Proceeding in the same way as done for the Coulomb problem, 
we see that for a point just below the branch cut, 9i ~ 27i and 6*2 ~ tt. 
Substituting these values in ()3.14p . one gets 



Pc 



A 



(3.15) 
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For pc to be positive just below the branch cut, we need to fix A = 1. Thus 
in the entire complex y plane, 



Pc = z^^==^(rir.)V^exp ( z^^^i^ ) , (3.16) 



which in turn gives 

Xc = -^^^(nr,)^/^expU^l^l (3.17) 



3.2.3 Residues at the fixed poles y = ±1 

At y = -1 : We calculate the residue of Xc at the fixed pole y = —1. In the 
complex plane the point y = ~1 corresponds to ~ ~ tt, which gives 



Xc = \_^,, {rir2r' - TTTT^. (3.18) 
The residue of Xc at ?/ = — 1 is 

(rir2)^/^ + -, (3.19) 

which is a positive quantity. Now, we evaluate the residue of x at ?/ = — l,by 
doing Laurent expansion of x around ?/ = — 1, as 

X = -^ + «o + ai(|/ + l) + (3.20) 

Substituting this in (|3.7|) and comparing the coefficients of j^j^^-, one obtains 
the quadratic equation 

, 'ia^n^-{2A-.2Bf-,A{A^B)an 
(6i) - Ti\ + ^^-^ = 0, (3.21) 
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which gives the two values of b'l as 

{A + Bl n (A + B) 3h 

2a + 4 ' 2a + 4 " ^^'^^^ 

Since the residue of Xc at ?/ = —1 is positive, ()3.1U|) . imphes that we select 
the value of b'l, which is positive in the limit h 0. Thus, we pick 

r ^ + f if {A + B)>Q 
b[ = (3.23) 
l-^ + f if iA + B)<0. 

From the above equation, we see that both the values of b'l are acceptable, 

but in two different ranges. Comparison of these ranges with the ranges 

given in (|3.2|) and (|3.3|) . show that they correspond to the different phases of 

SUSY. 

At y = 1 : 

Repeating the above process at the other fixed pole, one obtains the residue 
of Xc at ?/ = 1 as 

y/E + A^ , ,1/9 ^ , , 
(rir2)^/^ + -, (3.24) 

which is positive. The two values of the residue of x are 

[A-B^ n {A-B) sn 

Comparing equations ()3.24|1 and ()3.25|1 . one sees that each value of bi is 
acceptable in different ranges, i.e 

^ + 1 if {A-B)>0 
bi = { (3.26) 
-^ + f if {A-B)<0 
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But, the present case of interest is {A — B) > 0, for which the choice of 
residue is 

(A-B) h , , 

The range {A — B) > corresponds to both phases of SUSY, as seen from 
fl3.2p and ()3.3|) . Thus, from the above discussion, one can see that by applying 
the boundary condition ()3.1()|1 properly, one can obtain both the ranges of the 
parameters and simultaneously select the values of the residues which give 
physically acceptable solutions. We now proceed to obtain the eigenvalues 
and eigenfunctions for both the phases. 

3.2.4 Case 1 : {A - B) > 0, {A + B) > 

For the above range of parameters, x iii ^3.8|1 becomes 

, , (A~B h\ 1 fA + B h\ 1 ^P'(y) ^ , 

where, the suitable values of b'l and bi are taken from ()3.23p and ()3.27p 
respectively. Substituting x from the above equation in ()3.7p and putting 
^ = 1, we compare the constant terms and the coefficients of l/y^ terms 
individually, in the limit of y going to infinity and obtain C = and the 
energy eigenvalue as 

En = {A + naf - A^. (3.29) 
The equation for Pn{y) assumes the form 

(1 - y^)P" + + l)y)P' + n(n + /i + iy)P^{y) = (3.30) 
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with 

A~B A+B 
//= , z/= , (3.31) 

a a 

is seen to be the Jacobi differential equation and hence Pn{y) coincides with 
p^-i/2,u-i/2(^y'^ and the bound state wave function for the phase in which 
SUSY is exact, is 

V^„(l/) =iV(l-y)'^/2(l + y)-/2pr^'^-^(^). (3.32) 

The energy eigenvalues in ()3.29p and the wave functions agree with the so- 
lutions given in pp. 

3.2.5 Case 2 : (A- S) > 0, (A + S) < 

The right choice of the residues for this range, as shown in §3. 2. 3 is, 

. _(i±i?) ^ |. ,3.33, 

2a 4 2a 4 

Substituting x from (j3.8|) . with the above values of the residues in (j3.7p with 
^ = 1, one obtains the energy eigenvalue as 

E=(^B-{n + ^)a^ - (3.34) 

and the differential equation for Pn{y) turns out to be 

{l-y^)P:{y) + {iu-fi + l)-yiu + fi + 2))P:,iy)+n{n + fi + u+l)P4y) = 

(3.35) 

and the bound state wave function is found to be 

My) = N{1 - i/)t(l +y)i+^Pr^'^"^(y) (3.36) 
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where 

A-B 1 A+B 

/i = , z/ = (3.37) 

a 2 a 

The Scarf - I potential can be related to the Poschl - Teller - I potential by 

redefinition of the potential parameters as 

x^x + ^, A = 6 + (3, B = 6-(3. (3.38) 

With this redefinition of the parameters, we see that the bound state wave 
functions in ()3.36|) match with the bound state wave functions of the broken 
SUSY phase of the Poschl - Teller - I potential given in [7]. 

Thus, we see that the QHJ formalism in one dimension gives the accu- 
rate expressions for the bound state wave functions, when there are differ- 
ent phases of SUSY. It may be remarked here that, in the range A — B < 
0, A + B < 0, SUSY is exact but the roles of if_ and iJ+ are interchanged. 
In the range A — B > 0, A -\- B < 0, SUSY is again broken. 

3.3 Phases of a periodic potential 

In this section, we study the Scarf potential 

Vi^) = -- (3-39) 
sm (^) 

with a as the potential period. The Schrodinger equation for this potential, 
with h = 1, is 



Chapter 3. Study of potentials exhibiting different spectra 



46 



with A and s, defined by 



2 2mEa?' 



, 1 2n 2mVQa^ 
- - s = — 



(3.41) 



This potential is known to exhibit two different types of spectra for different 
ranges of s, unlike the Scarf - I potential, which exhibits different bound state 
spectra for different ranges. In the range s > 1/2, one can see from the plot 
in fig. 3.2 



with s = 2 and a =1, allows ordinary bound stat 
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that the potential is an array of infinite potential wells. Hence, the particle is 
confined to one well and the wave function should vanish at a; = ±a. Thus, 
in the above range, the potential has bound state spectrum. Prom the above 
plot, one can also see that the minima of the potential is at zero and hence 
the energy eigenvalues are always greater than zero. For s, in the range 
< s < 1/2, one can see from the plot in fig 3.3, 
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that the potential is similar to that of the potential in a crystal lattice and 
hence there are bands of allowed ranges. The particle can escape to infinity 
and the wave function need not vanish at x = ±a, but it should not become 
infinity anywhere. Unlike the previous case one cannot make any statement 
on the energy eigenvalues just by looking at the plot of the potential. But 
in the following sections, we will see that the energy is greater than zero. 

For this potential, we will show that one can obtain both the ranges and 
their corresponding solutions, by doing a QHJ analysis. The corresponding 
QHJ equation for the Scarf potential with p = —iq is 



vr 



+ q' + 



Performing the change of variable 



sin'(^) 



0. 



y = cot 



TTX 



(3.42) 



(3.43) 



we make use of the results in §2.4. The expression for df/dx, expressed in y, 
becomes F{y) = + y'^)/a, which when substituted in ()2.70p gives the 
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transformation equations, 

g = --(l + y^)0; = x--^. (3.44) 

Thus, the QHJ equation ()2.71|) for x, becomes 

3.3.1 Form of the QMF x 

The QMF has n moving poles with residue one on the real line, corresponding 
to the nodes of the n*'^ excited state wave function. From ()3.45p . one can see 
that X has fixed poles at y ±i. Hence, one can write x the rational form 

as 

,.J^ + JL + 5^ + C. (3.46) 

y-t y + t Pn[y) 

where hi and h\ are the residues at y = i and y = —i respectively and Pn{y) is 
an n*^ degree polynomial. C represents the analytic part of x and from ()3.45j] 
one can see that x is bounded for large y. Thus, from Louville's theorem, C 
is a constant and x is a rational function of y. The residues at the fixed poles 
y = ±i, when calculated in the same way as in the previous chapter gives, 

h = ; bi = . (3.47) 

We assume that x has finite number of singularities in the complex plane, 
which is equivalent to saying that the point at infinity is an isolated singu- 
larity. As mentioned in the end of §2. 4, one can make use of the fact that 
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for a rational function, the sum of all the residues is equal to zero, to find 
the eigenvalues, which gives 

bi + b[+n = di, (3.48) 

where di is the residue at infinity and is calculated by taking Laurent expan- 
sion of X around the point at infinity as 

X = do + - + -,+ , (3.49) 

y 

which when substituted in the QHJ ()3.45|) . gives 

dj-di + i^- s^) = (3.50) 

from which, the values of di turn out to be 

d, = (3.51) 

Substituting the values of the residues in ()3.48|) . one obtains 

n = -^±sTX, (3.52) 

which gives the degree of the polynomial Pn{y) in ()3.46|) . From ()3.41|) . one 
can see that if E < 0, X becomes imaginary. Hence, (j3.52p is not satisfied. 
Thus, from the above equation, we have the condition E > 0, which in turn 
implies A > and real. Hence, one can say that for any range of s, the energy 
eigenvalues are greater than zero. With this condition on A, we now proceed 
to select the values of residues at the fixed poles and at infinity, which will 
give us the physically acceptable results. 



Chapter 3. Study of potentials exhibiting different spectra 50 
3.3.2 Choice of the residues 

First, we shall fix the value of the residue at infinity. From the prior discussion 
of the potential, we know that the wave functions should not become infinite 
anywhere, in particular, for x = ±a, i.e. the wave function tp{y) should be 
finite for large y. From ()2.74p . the wave function in terms of x becomes 

ipiy) = exp - Y^^) ■ "^^'^^^ 

For large y, the leading behaviour of x is obtained as % ~ ^ from ()3.49|) . 
which when substituted in fl3.53|) . gives 



i^iy) - [j [j - YT^) dy) (3-54) 



yd, 



(3.55) 



(y2 + 1)1/2 • 

Substituting the value of di from ()3.51|1 in the above equation, one obtains 

^(^) ~ (^2^,^ 1)1/2 - (3-56) 

For < s < |, from the above equation, one can see that ip ^ for both 
values of di and this range corresponds to the case, where the potential 
exhibits band structure. 

For s > 1/2, ip{y) 0, only if di takes the value 1/2 — s. Thus, we 
fix the value of cii = 1/2 — s. In this way, the two different ranges of the 
potential parameter s emerge simultaneously fixing the values of di. To select 
the values of bi and b[, we note that the bound state and band edge wave 
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functions of one dimensional potentials are non- degenerate, which implies, 
the wave functions are real. Hence, x must be real, for which we must have 

61 = b[. (3.57) 

With the above constraint, ()3.48|) becomes 

2bi+n = di. (3.58) 

Thus, we see that the property, the wave function should be finite for x —>■ 00, 

gives the values of di in the two ranges as follows 

r for < s < 1/2 

di = I (3.59) 
[ for s > 1/2. 

From the parity constraint, one obtains the restriction on the values of the 

residues at the fixed poles as 61 = b[. Using these results, we proceed to 

calculate the solutions for the two ranges. 

3.3.3 Case 1 : Band spectrum 

In the range < s < 1/2, we have seen that di can take both the values 
of the residues. Hence, taking all the possible combinations of the residues, 
keeping 61 = b[ and substituting them in ()3.48j) . we evaluate for n, which 
gives the degree of the polynomial Pn{y) and takes the integer values greater 
than or equal to zero. Thus, there will be four combinations forming four 
different sets with the expression for n, given in the fifth column of table 
3.1. Since n needs to be positive, from table 3.1, we pick only those sets 
which give a positive integer value for n. As seen in §3.3.1, A is a positive 
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Table 3.1: All possible combinations of residues for the range < s < 1/2 



Set 


bi 


b[ 


rfl 


n = di — hi — b[ 


remark 


1 


l-A 
2 


l-A 
2 




= 1/2-s 




A > (s + i) 


2 


l-A 
2 


l-A 
2 


di 


= 1/2 + s 


A + s-i 


A>-(s-i) 


3 


1+A 
2 


1+A 
2 


di 


= 1/2 - s 


-A-.-i 


not valid 


4 


1+A 
2 


1+A 
2 


di 


= 1/2 + s 


-A+s-i 


not valid 



real constant. Thus, only the sets 1 and 2 will give positive values for n and 
hence, the other two sets namely 3 and 4 are ruled out. So, taking the values 
of bi,b[ and di from the sets 1 and 2 and substituting them in ()3.58|) . and 
using the definition of A and s from (j3.41|) . we obtain the expressions for the 
energy eigenvalues corresponding to the two band edges of the n*^ band, as 

^2/1 



" 2ma2 



n -\ — ± s 
2 



(3.60) 



Here, correspond to the upper and lower band energies of the n*^ band 
and they match with the solutions given in |5 13] . 

In order to obtain the expression for the wave function, we substitute x 
from ()3.46|) in ()3.53|) with bi = b'l, which gives 



Pn{y) 



(3.61) 



To evaluate the expression for the polynomial, we substitute x from ()3.46|) . 
with hi = b'l, in the QHJ equation ()3.45|) . which gives the second order 
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differential equation 

Pniy) + [-^) Pniy) + 

1/4 - 4%^ + 26i(l-^^) + A^-l 

2/2 + 1 + (1/2 + 1)2 

Tlie sets 1 and 2 liave 6i = (1 — A)/2, whicli on substitution in tlie above 
equation, gives 

PM + (^1^) ^;(!') + '~lf\'l)'' P'-'y^ = »■ (3,63) 

From ()3.4ip and ()3.60|) . one can see that A lias two values X = n ± s + 
|. Substituting these in the above equation, one obtains two differential 
equations corresponding to the two energy eigenvalues as 

{y' + l)P:{y) + (1 - 2n + 2s)yPl,{y) + n{n ± 2s)P„(y) = 0. (3.64) 

Defining y = it, the above equation takes the form of the well known Jacobi 
differential equation 

(l-t2)P;;'(t) + (z/i-z/2-t(z/i+z/2+2))P^(t)+n(n+z/i+i/2 + l)P„(t) = 0, (3.65) 

with Ui = 1/2 = —n + s — 1/2, for the corresponding two A values. Hence, 
the polynomial P„(t) will be proportional to the Jacobi Polynomial and the 
expression for the two band edge wave functions for the n^^ band will be 

^P{y) = {y-' + lf^P:^'-\-iy), (3.66) 

with respective z/i, z/2 values corresponding toA = n±2 + l/2. 



Pn{y) = 0. (3.62) 
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3.3.4 Case 2 : Bound state spectrum 

We proceed in the same way as done in case 1 i.e. take all possible combi- 
nations of bi,b[ and di, keeping bi = b[ in ()3.58|) . Since, di can take only 
single value 1/2 — s, only two sets are possible. Out of these, the set corre- 
sponding to 6i = 6'^ = (1 — A)/2 alone, will give a positive value for n. Thus, 
substituting these values of residues in ()3.58p . one obtains the expression for 
the energy eigenvalue as 

where n can take values 0, 1, 2... and is left with the Jacobi differential 
equation in terms of t as 

(1 - t')P:{t) - 2t{-n - s + ^)P^(t) - n(n + 2s)Pn{t). (3.68) 

The expression for the wave function from ()3.53|1 becomes 

^{y) = {y'+l)-^P^^^-^^{-ty), (3.69) 

where Si = S2 = —n — s — 1/2. 

Note that by subtracting off the ground state energy and redefining a = 
— 1, A = —{s + 1/2) and x ^ vrx/a in 

V = A{A - a)cosec ^ax - A^, (3.70) 

the supersymmetic potential P, one can see that the above expression and 
its solutions match with those obtained here. 
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Chapter 4 



Exactly Solvable Periodic Potentials 
4.1 Introduction 

This chapter is devoted to the study of potentials, which are ES and peri- 
odic in nature. It is well known that periodic potentials play a significant 
role in condensed matter physics. The energy spectrum of periodic poten- 
tials is unique due to the existence of energy bands and the solutions of the 
Schrodinger equation have the Bloch form: 

ip{x) = u{x) exp{ik.x) (4.1) 

Here, u{x) has the periodicity of the potential function. These features of the 
periodic potentials are usually illustrated using the Kronig - Penney model 
PP, in which the crystal lattice is approximated by a periodic array of square 
wells. Imposing the Bloch condition on the solutions of the Kronig-Penney 
model, one is led to a transcendental equation, which has to be solved, in 
order to obtain the band edge solutions. 

Recently, there has been considerable interest in the study |2 - [7] and 
construction [U |H1 E] of new periodic potentials, especially those belonging 

56 
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to the elliptic class. The manifestation of periodic lattices in BEC has also 
renewed interest in this area |l()pilj . This gave us enough motivation to check, 
if we could apply the QHJ formalism to these potentials. In the previous 
chapter we have discussed the Scarf potential [121112], V{x) = ~ 2msin^(^) ' 
which exhibits band structure in a finite domain < s < 1/2 and obtained 
the band edge solutions of this potential. 

Here, we study the ES case of Lame and the associated Lame poten- 
tials, which belong to the class of elliptic potentials and obtain the band 
edge solutions ^1]. These potentials are usually expressed in terms of the 
Jacobi elliptic functions sn (x, m), cn (x, m) and dn (x, m) [13 CHj, where the 
parameter m is known as the elliptic modulus, whose value lies between 
and 1. 

The Lame potential 

y(x) = j(j + l)msn^(x, m), (4.2) 

with j being an integer, has 2j + 1 bands followed by the continuum. This 
potential has been extensively studied in It is interesting to note 

that the Schrodinger equation with the Lame potential, 

d tlj{x) ^ 2m ^ ^ l)msn^(x,m))V^(x) = (4.3) 

is a result of separating the Laplace equation in the ellipsoidal coordinates 
and is referred to as the Lame equation [T7[ I19j. This has unique analytical 
properties fTj, which render it useful in the study of topics ranging from 
astrophysics [201 -EH], condensed matter physics [2Z1 to chaos |28|29j . 
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In comparison to the Lame potential, the associated Lame potential 



V{x)=pmsn [x,m) + qm- — 57 r, (4.4j 

dn (x, m) 

with p = a{a + 1) and q = b{b + 1), did not have a systematic study till 
recently and has a few scattered references in the mathematical literature. 
This potential is ES, when a = b = j with j being an integer, and QES when 
a ^ b with a, b being real. A systematic study of the band structure of this 
potential has been done for various cases of a and b by Khare and Sukhatme 
PU] and the algebraic aspects of the QES associated Lame potential have 
been investigated in ^11- 1^3]. 

The band edge eigenf unctions and the eigenvalues of both these potentials 
have been found using the conventional techniques for solving the ordinary 
differential equations f7[E]. One can obtain the solutions, using the same 
method used to solve the band structure case of the Scarf potential in the 
previous chapter. 

In §4.2, we show, how one can obtain the form of the wave functions for 
the Lame potential for j being an integer, followed by §4.3, where we take 
up a special case of j = 2 and obtain explicit band edge wave functions and 
energies. In §4.4, the form of the wave functions for the general associated 
Lame potential (ES case) is derived, followed by the specific example, where 
of a = 6 = 1, in §4.5. 
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4.2 General Lame potential 

The QHJ equation in terms of q for the Lame potential is 

+ ^ + E + l)msn^(x,m) = 0. (4.5) 

As discussed in §2.4, we perform a change of variable 

t = sn (x, m), (4.6) 



which gives F{t) = J {1 — t^)(l — mt^). The elliptic function is periodic with 



sn (x, m) = sn (x + 4K{m),m), where K{m) is dtj y (1 — t^^{X — vrd?\ the 

complete elliptic integral of the first kind [El CEl inHj- We would like to 

point out here, that the Lame equation can be written in five forms, namely, 

two algebraic, one trigonometric, one Weierstrassian and one Jacobi ^7] 

depending on the change of variable. Of all these change of variables, we 

found ()4.6p . which gives the Jacobi form of this Lame equation, to be the 

best suited for our method. It enables one to write the QHJ equation in 

a form, which can be easily analyzed using the QHJ formalism. Another 

added advantage of this particular choice is, that it maps half of the period 

parallelogram < x < IKim) of the Jacobi elliptic sn (x, m) function to the 

upper half complex plane jHS]- Using the transformation equations ()2.70|1 . 

one gets the equation for x as 

2 (mt)^ + 2m t' + 2 2E - 2j(j + l)mt^ - mt=^ _ , . 

^ rft 4(1- mt2)2 + 4(1-^2)2+ 2(l-t2)(i_mt2) ^ K 

with 
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The following properties of the Jacobi elliptic functions [E] 

— sn (x, m) = cn(x, m)(in(x, m) (4-9) 
ax 

and 

sn ^(x, m) + cn^(a:, m) = 1, sn ^(x, m) + m(in^(a:, m) = 1 (4-10) 

were used. For all our further calculations we will regard x &s the QMF 
instead of p. Using (j4.7|) . in place of the original QHJ equation (j4.5p . one 
can obtain the expressions for the wave functions by analyzing the singularity 
structure of x- 

4.2.1 Form of QMF x 

Equation ()4.7|1 shows that x has fixed poles at t = ±1 and t = ±l/-ym. 
We assume, there are a finite number of moving poles in the complex plane. 
Hence, we make an assumption that the point at oo is an isolated singular 
point and there are no singularities of the QMF except for those mentioned 
above. Therefore, proceeding in the same way as done in previous chapters, 
we can write separating it into its singular and analytical parts as 

Here, Q{t) is the analytic part of x ^^id the rational terms represent the 
singular parts. The coefficients bi, b[ and di, d[ are the residues at t = ±1 
and if: = ±l/^/rn respectively. The term P^(t) / Pn(t) represents the singular 
part coming from the n moving poles with residue one. Q{t) represents the 
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analytic part of x and from ()4.7j) . one can see that x is bounded at infinity. 
Hence, from Louville's theorem, Q{t) is a constant. 

The residues 61, b[, di and d[ at the fixed poles can be determined by 
substituting the Laurent expansion of x around these poles in ()4.7j) . The two 
values of the residues at all the fixed poles turn out to be 



-, -. (4.12) 
4' 4 ^ ^ 



Note that all the residues are independent of the potential parameter j. 
Hence, irrespective of the value of j in the potential, the residues at the fixed 
poles will take only the above values for any Lame potential. 

We need to select the values of the residues, which give rise to acceptable 
wave functions. Similar to the case of Scarf potential, there is no way of 
selecting the residues, hence, we accept both the values i.e. 3/4 and 1/4 for 
each of the residues in the present case . The only requirement that restricts 
possible combinations of values is that of parity, which arises due the fact 
that the band edge wave functions are non - degenerate, iplt) having definite 
parity implies ip{—t) = ±'?/'(t), which in turn gives x{~t) = ~xit)- This 
condition rules out the possibilities bi 7^ b[ and di 7^ d[. So, the possible 
combinations of 61, b[, di and d'l, which are allowed, are those which satisfy 
bi = b[ and di = d'^. 

We have assumed that x has finite number of singularities in the complex 
plane and hence, the point at infinity is an isolated singularity of x- The form 
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of ()4.7|1 suggests that x is bounded at oo. Therefore, x can be expressed as 

x(t) = Ao + y + ^ + .... (4.13) 

and the coefficients A^'s are fixed using the QHJ equation ()4.7j) and one gets 
Ao = , and two values for Ai as 

Ai=j + 1, -J. (4.14) 

Making use of the fact, that for a rational function, the sum of the residues 
is equal to zero, we get 

2bi + 2di + n = Xi, (4.15) 

where, bi = b[ and di = d[ have been used. As the left hand side is positive, it 
is clear that the only choice Ai = j + 1 is consistent with the above equation, 
thus fixing the value of residue at infinity. 

Thus from ()4.15p . one gets the expression for n as 

n= j + l-2bi-2di, (4.16) 

which gives the degree of the polynomial Pn{t). So making use of all the 
residue values at the fixed poles keeping the constraint bi = b[ and di = d[, 
one obtains four different sets and corresponding four different values of n 
from ()4.16|1 . as listed out in table 4.1 



Table 4.1: Number of zeros of the wave function corresponding to the number of moving poles n of x, form 
of ip{x), number of linearly independent solutions and the total number of zeros of ip{x). 



9 



set 


h 


di 




n 


ijj in terms of x 


Number of LI solutions 


Total zeros 














J = 27V + 1 


j = 2N 


of ip{xy 


1 


1/4 


1/4 




j 


Pj (snx) 


N + 1 


N+1 


j 


2 


3/4 


1/4 


j 


- 1 


cn xPj_i{snx) 


N+1 


N 


j 


3 


1/4 


3/4 


J 


- 1 


dn xPj_i(snx) 


N + 1 


N 


J-1 


4 


3/4 


3/4 


j 


- 2 


cn xdn xPj-2 (sn x) 


N 


N 


J-1 



EC 
O 

CO 



§ 

Co 



These constitute of both real and complex zeros. 



CO 
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4.2.2 Explicit forms of the wave functions 

One can obtain the form of the wave functions using ()2.74|) . Substituting 
the relations given in ()4.8j) in ()2.74j) . one gets the expression for the wave 
function in terms of x as 

*W=exp(/(,w4(^ + ^))*). (4.17) 

Substituting x from ()4.1H) gives 

which when simplified and written in terms of the original variable x, gives 

)°(dnx)^P„(sna;), (4.19) 

where, a = ^^^^ ~ ^^2~^" -^^^ ^ake of simplicity, the elliptic 

modulus m of the Jacobi elliptic functions will be suppressed here onwards. 
The four different combinations of the residues (given as sets 1 to 4 ) give 
rise to four different forms of the band edge wave functions as listed in the 
table 4.1. 

The parity constraint x(~^) = ~x{^) restricts the polynomial P„(t) to 
have either only odd or only even powers of t = snx. In the sixth and the 
seventh columns of table 4.1, the number of linearly independent solutions, 
for the two cases, j being odd and j being even, are given in terms of a 
positive integer N. For odd j, = (j — l)/2 and for even j , N = j /2. In 
both the cases, the total number of solutions for a particular j is equal to 
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2j + 1. It is easy to see that the four forms and the number of solutions of 
a particular form obtained here are in agreement with those already known 
[T7j . For a given set of bi and di, n is fixed using ()4.16|) and the differential 
equation for the unknown polynomial Pn(t) can be obtained by substituting 
x(t) from ()4.1H) in (jUTj), which gives Q = and 



For each set of residues given in table 4.1, this differential equation is equiv- 
alent to a system of n linear equations for the coefficients of different powers 
of t in Pn{t). The energy eigenvalues are obtained by setting the correspond- 
ing determinant equal to zero. We illustrate this process by obtaining the 
eigenvalues and eigenfunctions explicitly for the case j = 2. 

4.3 Lame potential with j = 2 

Using the procedure described in the previous section, we obtain the eigen- 
values and the eigenfunctions for the supersymmetric potential [H] 




(4.20) 



where 



Git) 



t\Abj - 2bi + 1/4) + 1/2 - 2bi 
(t2 - 1)2 

{mtfiMj - 2di + 1/4) + m/2 - 2mdi 
^ (mt2 - 1)2 

2E - 2j{j + l)mt2 + (IQmbidi - l)mt^ 
2(t2 - l)(mt2 - 1) 



V- (x) = 6m sn — 2m — 2 + 25, 



(4.21) 
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where 6 = yl — m + m?. This potential is same as the Lame potential in 
(14. 2j) with j = 2, except for an additive constant 26 — 2m — 2, which has been 
added to make the lowest energy equal to zero. The equation for x(t) is 

^ dx {mtf + 2m t^ + 2 2E + 4{m + 1 - 6) - 13mt" _ 
^ rft 4(l-mt2)2 + 4(l-t2)2+ 2(1 -t2)(i_,„t2) u. ^ . J 

From ()4.16|) . the number of moving poles n, for j = 2, are 

n = 3 - 26i - 2di, (4.23) 

where the values of bi and di are obtained from and ()4.12|) . For each set of 
bi, di and n, one can write the form and the number of linearly independent 
solutions by substituting j = 2, which gives = 1 in table 4.1. Hence for 
the case j = 2, the form and the number of solutions are as given in table 
4.2. 

To get the unknown polynomial part Pn{t) of the wave function and the 
band edge energies, one needs to substitute the different sets of bi, di and n 
from table 4.2 in the differential equation 

P:{t) + 4t + PLit) + G{t)Pn{t) = 0, (4.24) 

where 

t\Abj - 26i + 1/4) + 1/2 - 2bi 



Git) 



+ 



(t2 - 1)2 

(mt)2(4rf2 - 2di + 1/4) + m/2 - 2mdi 
(mt2 - 1)2 

2E + 4(m + 1 - 5) + (166irfi - 13)mt'^ 
2(t2-l)(mt2-l) ■ 
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Table 4.2: The form of the wave functions and the number of hnearly inde- 
pendent solutions for the Lame potential with j = 2. 



set 


bi 


di 


n 


ip in terms of x 


Number of LI solutions 


1 


1/4 


1/4 


2 


P2(sna;) 


2 


2 


3/4 


1/4 


1 


cn xPi (sn x) 


1 


3 


1/4 


3/4 


1 


dn xPi(snx) 


1 


4 


3/4 


3/4 





cnxdn xP(){snx) 


1 



Thus, for various sets of residues, one has the band edge energies and the 
wave functions as: 
Set 1 : 6i = 1/4 , rfi = 1/4, n = 2 
Taking P2 = At^ + Bt + C, the parity constraint implies 



B = 0. 



(4.25) 



Substituting bi, di and P2 in ()4.24|) . one gets a 2 x 2 matrix equation for A 
and C as follows 



E -2m -2 -26 



-6m 



A 



E + 2m + 2-26 / I C 



0. 



(4.26) 



Equating the determinant of the matrix in ()4.26|) to zero, one gets the two 
values for the energy as 



Ei = 4:S, Eii = 0, 



(4.27) 
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which in turn give the polynomial as 

Pj{t) = m + 1 - 5 - 3mt^ , Pji{t) = m + 1 + 5 - 3mt'^ . (4.28) 
From ()4.19|) . we see that the band edge wave functions in terms of x will be 
ipi{x) = m + 1 — 5 — 3m sn , ^ii{x) = m + 1 + 5 ~ 3m sn ^(a;). (4.29) 
Set 2 : 6i = 3/4 , rfi = 1/4, n = l 

Taking Pi = t — to and substituting in ()4.24p , one gets to = and the band 
edge eigenvalue as 

E = 2S -m + 2. (4.30) 
The band edge wave function becomes 

ip{x) = cnxsnx. (4-31) 

Set 3 : 6i = 1/4 , rfi = 3/4, n = l 

Proceeding in the same way as in set 2, one gets 

E = 25 + 2m-l, iIj{x) = dnxsnx. (4.32) 

Set 4 : 6i = 3/4 , rfi = 3/4, n = 

In this case, Po(t) is a constant and one gets the band edge energy and the 
wave function to be 

E = 2S — m — 1 , il){x) = cnx dnx. (4.33) 

Thus, one obtains five band edge wave functions and their corresponding 
energies, which agree with those given in [S]. 
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4.4 Associated Lame potential 

The associated Lame potential is ES when a = b = j. In this case, the 
potential expression becomes 

(cn \ 
sn^x + —^]. (4.34) 
dn X y 

The QHJ equation with the associated Lame potential is 

q^ + q' + E- mj(j + 1) (sn^x + ^] = 0. (4.35) 

\ dn a; / 

Doing the change of variable t = snx and proceeding in the same way as in 
the case of general Lame potential described in §4.2, one gets the equation 

for x{t) as, 

2 dx (mt)^ + 2m(l-2j(j + l)) 2 + 
^ dt 4(1 - mt^y 4(1 - ^2)2 

+^^^^^!i^±M^.O. (4.36) 
4.4.1 Form of QMF (x) 

Similar to the Lame potential, x(t) has fixed poles at t = ±1 and t — ±l/-\/m 
along with n moving poles in the entire complex t plane. Further, we assume 
that the point at infinity is an isolated singularity, x has no other singular 
points, except at t = ±1 and t = ±l/y^. Following the same procedure 
used in previous chapters, one obtains the residues at i = ±1 as 



Chapter 4. Exactly Solvable Periodic Potentials 70 

which is independent of j. The residues at t = ±l/^/m turn out to be 

3 + 2j l-2j 3 + 2j l-2j 

di = , d, = (4.38) 

which are j dependent. Knowing the singularity structure of x, one can write 
it as 

y- + ^'^ I d[ P'St) 

vahd for all t similar to the Lame potential. Assuming that the point at infin- 
ity is an isolated singularity, we calculate the residue at infinity by expanding 
X in Laurent expansion for large t as 

X = Ao + ^ + ^ + .... (4.40) 

and substituting it in the the QHJ equation (|4.36p . one obtains the two values 
for Ai as 

Ai = J + 1, -J. (4.41) 

For a rational function, the sum of all its residues is equal to zero. Hence, 
for X 

hi + h\ + di + d[+n = \i. (4.42) 

As in the Lame case, the parity requirement implies hi = 6'^ and di = d[ and 
we get 

26i + 2c/i + n = Ai. (4.43) 

In this case, both the values of Ai are allowed because di and d[ can take 
negative values unlike the case of Lame potential, where Ai = —j was ruled 



Chapter 4. Exactly Solvable Periodic Potentials 



71 



out. Thus, one has two cases Ai = j + 1 and Ai = —j. Note that the 
potential is invariant under the transformation j —>■ — j — 1 and one has 
two different values of j, one positive (j = j' > 0) and another negative 
(j = — / — 1 < 0), leading to the same expression for the potential described 
by ()4.4|) and hence, the same answers for the eigenf unctions and eigenvalues. 
Therefore, it is sufficient to restrict j to positive values alone. We consider 
only the case Ai = j + 1 and neglect the other value Ai = — j. For \i = j + 1, 
the values of n, which describe the number of zeros, for different combinations 
of bi and di values, are given in table 4.3. 

The sets 3 and 4 give n = —1 and n = —2 respectively and will not be 
considered, as n should be greater than or equal to 0. The sets 1 and 2 will 
give positive values of n, only if j is positive. Since the singularity structure 
of the associated Lame potential is same as the Lame potential, the form 
of the wave functions will be same as that given in ()4.19p . As in the Lame 
case, the number of linearly independent solutions are given in the table 4.3 
in terms of N, for j being odd or even integer. 

4.5 Associated Lame potential with j = 1 

We perform the calculation with the supersymmetric potential 




dn X 



(4.44) 



corresponding to j = 1. 



Table 4.3: Values of n and the form of the wave functions, for different combinations of bi and di, for 
Xi=j + 1. 



set 


bi 


di 


n 


■?/'(a;) in terms of x 


Number of LI solutions 


Total zeros 












j = 2N + 1 


j^2N 


of ijj{x)'^ 


1 


1/4 


(1 - 2j)/4 


2j 


{dnx)~^ P2j{snx) 


2N + 2 


2N + 1 


2j 


2 


3/4 


(1 - 2j)/4 


2j-l 


cn a;(dn x)~^ P2j-i{sii x) 


2N + 1 


2N 


2j 


3 


1/4 


(3 + 2i)/4 


-1 










4 


3/4 


(3 + 2j)/4 


-2 











^These constitute of real and complex zeros. 



Chapter 4. Exactly Solvable Periodic Potentials 



73 



Proceeding in the same way as in the previous sections, the equation for 
x(t) is found to be 

2 dx {mty — Qm + 2 

^ ^ 4:{mt^ - 1)2 + 4(t2 - 1)2 

2E + A + 2m- 4V1 - m - hmt^ 

+ 2(i-,.e-)ii-P) = « (^-^i 

and the form of x is 

bi h\ di d'-. P'it) ^ 

y = — — H — H H + + C. (4.46) 

^ t-1 t+1 t-l/v/^ t + l/v/^ P„(t) ^ ' 



Substituting ()4.4fi|l in ()4.45|1 . one obtains C = and is left with the differen- 
tial equation 

Pnit) + 4t + ^^^^ K(t) + G(t)P„(t) = 0, (4.47) 

where 

t2(462 - 2bi + 1/4) + 1/4 - 2hi 



Git) 



+ 
+ 



(t2 - 1)2 

(mt)2(4rf2 _ 2di + 1/4) + -3m/2 - 2m(ii 
(mt2 - 1)2 

2E + A + 2m- Ay/l-m + (IQbidi - 5)mt^ 



2(t2-l)(mt2-l) 

Substituting j = 1 in table 4.3, one gets the band edge eigenfunctions 
and eigenvalues from set 1 and set 2 as follows. 
Set 1 : 6i = 3/4 , rfi = -1/4, n = 1 
This combination gives only one solution with 



(4.48) 
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Set 2 : 6i = 1/4 , = -1/4, n = 2 

This combination gives two solutions. The band edge energies are 

Ei = 0, En = 4Vl -m (4.49) 

and the corresponding wave functions are 

, / N 1/^1 a/1 — ^\ , / N 1 /i Vl — m\ ,^ 

^/'/(a;) = — dnx H , ipn[x) = — dnx , (4.50) 

m \ anx J m \ dnx/ 

which match with the wave functions in [S]. Thus, in this chapter, we have 
shown that the QHJ formahsm successfully yields the band edge eigenvalues 
and eigenfunctions for ES periodic potentials. 

For all the models studied, the parameter n, which appears in the exact 
quantization condition ()2.4()j) described in chapter II, gives the number of 
moving poles of the QMF, which in turn correspond to the zeros of the wave 
function. For ES models [HZj, these moving poles were found to be confined 
to the classical region alone. Thus, the n^^ excited state has n real zeros. 
The parameter n characterizes the energy level and increasing n increases 
the energy. 

For the QES models, the same n appears as a parameter in the potential 
and of the possible states, only a few states, which are determined by n, can 
be obtained analytically. The QMF for all these states will have the same 
number of moving poles, which are both real and complex. Correspondingly 
the wavefunctions for each state will have the same number (= n) of zeroes 
(both real and complex), of which the number of real zeros is in accordance 
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with the oscillation theorem. The number of real zeros increases as one goes 
from the lowest state to the highest state possible. 

On the other hand, the exactly solvable periodic potential models show 
quite a different kind of distribution of the moving poles in the complex 
plane. Like QES models, the moving poles of the QMF consist of both real 
and complex poles. But unhke the QES models, where the number of poles 
of QMF for all the known states is same, for ES periodic potential models one 
finds groups of solutions with number of moving poles remaining same within 
a group, but varying from one group to another group of solutions. This point 
becomes clear from the form of the solutions listed in table 4.1 for the Lame 
potential. Different groups of solutions are precisely the different sets of the 
solutions hsted in table 4.1. In each set, the degree of the polynomial Pn(sn x) 
is same and the total number of zeros of the band edge wave functions in 
the interval < x < 2K{m) are same. For example, for the + 1 linearly 
independent solutions belonging to set 1, the total number of zeros of the 
wavefunctions will be j. However, for the solutions belonging to set 2, there 
will be J — 1 zeros of P„(sn x) and one zero corresponding to cn thus a total 
of j zeros. For all the sets, the total number of zeros are given in the last 
column of table 4.1. However, the number of real zeros in a group increases 
with the increasing energy and correspondingly the number of complex zeros 
will decrease. For all values of j, we will have utmost four such groups. A 
similar pattern is observed for ES associated Lame potential. 
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Chapter 5 



Periodic Quasi Exactly Solvable Models 
5.1 Introduction 

In chapter I, we have mentioned that the QHJ analysis of ordinary QES mod- 
els i.e. potentials with ordinary bound state spectrum, successfully yielded 
the QES conditions along with the eigenvalues and the eigenf unctions |ll2j . 
The present chapter is a study of QES periodic potentials. Here, we discuss 
the QES case of the associated Lame potential [Hj. The QHJ analysis works 
on certain assumptions on the singularity structure of the QMF, hence, it 
is essential to see the validity of these assumptions by analyzing as many 
different models as possible. Therfore, we study the QES associated Lame 
potential, even though the method used is very much similar to that used in 
the previous chapter. 

The expression for the associated Lame potential is 

V [x) = a[a + l)msn (a;, m) + 0(0 + Ijm- — 27 r- (5-1) 

dn (x, m) 

Physically, this potential is described as a periodic lattice of period K{m) 
with the basis composed of two different atoms, which are alternately placed. 

80 
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For a systematic study of the associated Lame potential, one is referred to |H 
IS]. This potential is ES, when a = b and has been analyzed in the previous 
chapter, QES properties when a ^ b. There are three different QES cases 
depending on whether a and b are integers or non - integers in j3] and are 
given below. 

Case (i) : For a, b being unequal integers with a > 6 > 0, there will be a 
bound bands followed by a continuum band, if a — 6 is odd(even) integer, 
then there will be b doubly degenerate band edges of period 2K{AK) with 
the corresponding band gaps zero, which cannot be obtained analytically. 
Case (ii) : For a, b being half - integers, but a + b and a — b being integers, 
there will be infinite number of bands with band edge wave functions of 
period 2K or 4K depending on whether a — 6 is an odd or even. Of these 
infinite bands, there will he a — b bands whose band edges are non-degenerate 
with period 2K{AK) and b + ^ doubly degenerate states of period 2K{4K) 
and can be obtained analytically. 

Case (iii) : For a being an integer and b being a half integer or vice versa, 
one can obtain some exact analytical results for mid - band states. 

For all these above mentioned cases, the information regarding the exis- 
tence of the missing states, mid band states etc. can be acquired by appealing 
to the oscillation theorem [HIIZI- 

We proceed with the analysis of the associated Lame potential with the 
same change of variable and the transformation equations used in chapter 
IV. With the same assumptions on the singularity structure of p, we proceed 
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to obtain the QES condition and the forms of the band edge wave functions 
for the general associated Lame potential in the next section. In section §5.3, 
we analyze case (i), where both a, b are integers taking the values 2 and 1 
respectively. Case (ii), where a and b are both half integers with the values 
7/2 and 1/2 respectively, is analyzed in section §5.4. 

5.2 QES condition and the forms of the wave functions 

The QHJ equation for the associated Lame potential, putting h = 2m = 1, 
is 

Ccn ^ fx) \ 
E - a{a + l)msn^{x) -b{b + l)m — • (5-2) 
dn [x) J 

Note that with the transformation b —b — 1 or a —a — 1, the potential 
does not change. Hence, for all our studies, without losing generality, we 
take a, b positive and a > b. Proceeding in the same way as in chapter IV, 
defining p = —iq and using the change of variable t = sn (x), one obtains the 
equation for x as 

2 dx mH'^ + 2m{l - 2b{b + 1)) 2 + t^ 2E - mt^ {I - 2a{a + I)) _ 
^ 4(l-mt2)2 +4(1-^2)2+ 2(l-t2)(l_r7it2) 

(5.3) 

with the help of the transformation equations 

Like all the earlier investigations, we shall treat x the QMF and ()5.3p as 
the QHJ equation. 
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5.2.1 Form of QMF (x) 

From ()5.3j) . we see that x has fixed poles att = ±1 and t = ±:^- In addition 
to the fixed poles, the QMF is assumed to have finite number of moving poles 
and no other singular points in complex plane. Hence, one writes x as a sum 
of the singular and analytical parts as follows 

where bi, b'l and di, d[ are the residues at t = ±1 and t = respectively, 
which need to be calculated. Pn(t) is an n*'^ degree polynomial with ^ = 
J2k=i being the summation of terms coming from the n moving poles 
with residue one. The function Q{t) is analytic and since x is bounded at 
infinity, from Louville's theorem, it is a constant, say C. The residues at the 
fixed poles can be calculated by taking the Laurent expression around each 
individual pole and substituting them in ()5.3|1 as in previous chapters. The 
two values of the residues bi ,b[ at t = ±1 are 

and at t = ±-^, one has the residue values 

3 6 16 ^, 3 b 1 b , , 

di = - + - , d\ = - + - , . 5.7) 

42'42 ^42'42 ^ ' 

As in chapter IV, we consider both the values of the residues. We demand 
bi = b[ and di = d\, a condition coming from the parity constraint x(t) = 
— We assume that the point at infinity is an isolated singularity and 
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hence using the fact, that the sum of all the residues of a rational function 
is zero and the restriction bi = b[ and di = d[, we obtain 

26i + 2di + n = Ai, (5.8) 

where Ai is residue at infinity. Since x has an isolated singular point at 
infinity, one can expand x Laurent series around the point at infinity as 

x(t)=Ao + y + ^ + .... (5.9) 

Substituting in (Q and comparing the various powers oft, one obtains 

Ai = a + 1, -a. (5.10) 

Taking various combinations of hi and di from ()5.6p and ()5.7|) . substituting 
them in (j5.8|) . one obtains the QES condition from each combination as given 
in table 5.1 for A = a+1. Thus, one sees that all the allowed combinations of 
residues give one of the forms of QES condition where n = 0, 1, 2.... Note 
that the other value of residue at infinity i.e. Ai = —a, when substituted 
instead of a + 1 in ()5.8|) . gives the QES condition for negative values of a, b 
i.e. for a —a — 1, b ^ —b — 1 in 6i, di. 

5.2.2 Forms of wave function 

Since the change of variable and the transformation equations are same as 
in chapter IV, we get the same expression for the wave function as in ()4.18p 
and ^^~T^ i.e. 

= exp (/ (x + \ (^ + ^)) dt) , (5.11) 
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Table 5.1: The quasi - exact solvability condition from the four permitted 
combinations of bi and di 



set 


bi 


di 


26i + 2(ii + n = Ai 


QES condition 


1 


3/4 


3,6 
i' 2 


2+b+n=a 


b — a = —n — 2 


2 


3/4 


1 b 
4 2 


1 — b + n = a 


a+b+l=n+2 


3 


1/4 


3 _i_ b 
i' 2 


1 + b + n = a 


b — a = —n — 1 


4 


1/4 


1 b 
4 2 


—b + n = a 


a + b = n 



which on substitution of x from ()5.5j) and written in terms of the original 
variable x, becomes 

)"(dnx)^P„(sna;), (5.12) 

where a = P = Hence, for each set of bi, di, one gets a wave 

function given by ()5.12j) . The degree n of this polynomial, which is obtained 
from ()5.8|) . as 

n = a + l-2bi~2di, (5.13) 

which is in terms of either a + b or a — b, as evident from table 5.1. The 
forms of the wave function can be found and are as given in table 5.2 and 
table 5.3, for the two different cases, when a + b and a — 6 are odd and even 
separately. 



Case I Both a + b, a — b are even : We introduce = ^ and M — 2=^, where M and N are integers, and 
obtain the forms of the wave functions in table 5.2, in terms of M and for the four sets of combinations 
of bi and di in table 5.1. 

Table 5.2: The form of the wave functions for the four sets of residue combinations when a + b and a — b 
are even and equal to 2N and 2M respectively. 



set 


bi 


di 


n = Ai - 26i - 2di 


n{M, N) 


wave function ip{x) 


LI solutions 


1 

2 
3 
4 


3/4 
3/4 
1/4 
1/4 


3,6 
4^2 

1 b 
4 2 

3,6 
4 2 

1 b 
4 2 


a-b-2 
a + b-1 
a-b-1 
a + b 


2M-2 
2N-1 
2M -1 
2N 


cn a; (dn x) ^^'^P2M-2{sn x) 
(f^^2iv-i(snx) 
(dn x)''~^^P2M-i{sn x) 

P2jv(sna:) 
(dna;)'' 


M 
N 
M 
N + 1 



Case II Both a + b, a - b odd : Introducing, N' ^ ^ and M' = where M' and N' are integers, we 
obtain the wave functions in table 5.3, in terms of M' and N' for the four sets of combinations of bi and di 
in table 5.1. 

Table 5.3: The form of the wave functions for the four sets of residue combinations when a + b and a — b 
are odd and equal to 2N' + 1 and 2M' + 1 respectively. 



set 


bi 


di 


n = Ai - 26i - 2di 


n{M, N) 


wave function il){x) 


LI solutions 


1 

2 

3 
4 


3/4 

3/4 
1/4 
1/4 


3 _|_ ft 

4 2 

1 h 
4 2 

3,6 
4 ~'~ 2 

1 b 
4 2 


a-b-2 
a + b-1 
a — b — 1 
a + b 


2M' - 1 

2N' 

2M' 
2N' + 1 


cn x(dn x)^"'"^P2M'-i (sn x) 
(dnxj^'+^PaM^sna;) 

-P2JV'+l(snx) 
(dna;)'' 


M' 
N' + 1 
M' + l 
AT' + l 
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From the forms of the wave functions in tables 5.2 and 5.3, one observes 
that the number of hnearly independent solutions is different for the two 
cases. The unknown polynomial in the wave function can be obtained by 
substituting x from ()5.5|1 in the QHJ equation ()5.3j) . which gives 

P^it) + 4P„(t) + + G(t)P„(t) = (5.14) 



where 



G{t) = 

+ 
+ 



t^(46j-26i + i)-26i + | 
(t2 - 1)2 

mH^{Adl - 2di + j) - 2mdi + ml 
(mt2 - 1)2 

2E + {Wbidi - 1 - 2a(a + l))mt'^ 
2(l-t2)(l_rnt2) • 



1-26(6+1) ■ 



The above differential equation is equivalent to a system of n linear equations 
for the coefficients of different powers of t in Pn{t). The energy eigenvalues 
are obtained by setting the corresponding determinant equal to zero. In the 
next section, we obtain the band edge wave functions for the associated Lame 
potential with a = 2, 6 = 1 

5.3 Associated Lame potential with a, b as integers 

For this case, we consider the associated Lame potential with a = 2, 6 = 1. 
Similar to the previous chapter, we work with SUSY potential 

2 

cn jC 

VI (x) = 6m sn ^x + 2m — ^ 4m. (5.15) 

dn X 
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This potential is same as (jS.lj) with a = 2 and 6=1, except that a constant 
has been added to make the lowest energy equal to zero. The QHJ equation 
in terms of x is 

2 , mH'^ - 6m 2 + ^2 2E + 8m - ISmt^ _ 
^ ^ ^ ^ 4(l-mt2)2 + 4(1-^2)2 + 2(l-t2)(i_,„t2) = 0- (5-16) 



Apart from n moving poles, x has poles at t = ±1 and t = il/i/m. As in 
the previous section, one can write x with the parity constraint as 

which gives the form of x iii the entire complex plane, where Pn{t) is yet to 
be determined. Note that for this potential the combination a + b and a — b 
are both odd i.e. 3 and 1 respectively. Hence, we use table 5.3 to obtain all 
the information regarding the residues at the fixed poles, number of moving 
poles of X) number of linearly independent solutions and their form for each 
set etc. by taking the values of a = 2, 6 = 1, which give M' = and A^' = 1. 
The unknown polynomial in the wave function can be obtained from ()5.14|1 . 
where G{t) for this potential turns out to be 

G{t) 



t2(462-26i + i)-26i + i 



(t2 _ i; 



2 



mH\4dj - 2di + i) - 2mdi - ^ 
^ (mt2 - 1)2 

2E + 8m + (166irfi - 13)mt2 
+ 2(l-t2)(i_^t2) • (5.18) 

Using ()5.14|) with G{t) from ()5.18|) and proceeding in the same way as is in 
the previous section, one gets the explicit expressions for the eigenf unctions 
and the eigenvalues as given in table 5.4. 



Table 5.4: The residues, the value of n, number of linearly independent solutions, the band edge eigenfunc- 
tions and eigenvalues are as follows. Here, a + b — 3 and a — b— 1, which give N' — 1 and M' — 0. 



set 


bi 


di 


n 


LI Solution 


Eigenfunction ipi^x) 


Eigenvalues 


1 


3/4 


5/4 


-1 








2 


3/4 


-1/4 


2 


2 


^(3msn - 2 ± V4 - 3m) 


5 - 3m ± 2 V4 - 3m 


3 


1/4 


5/4 





1 


dn^x 





4 


1/4 


-1/4 


2 


2 


^{Smsn'^x - 2 - m ± V4 - 5m + m^ 


5 - 2m ± 2v'm2 - 5m + 4 
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From the table 5. 4, we see that the first set of residues gives n = — 1, 
which will not be considered as n cannot be negative. Thus, this particular 
case of associated Lame potential has 5 band edge solutions, which can be 
obtained analytically out of an infinite number of possible states. 

5.4 Associated Lame potential with a and b as half 
integers 

The potential studied here is the supersymmetric associated Lame potential 

with a = 7/2,6 = 1/2, whose expression is 

T. 63 2 3 cn^a; ^ 29 ^ , 

V_ = — msn X -\ — m — ^ 2 m + Oq 5.19 

4 4 dn^x 4 

where = a/4 — 4m + 25m^. 

, , 2 + mH"^ -m AE + 8 + 29m - 450 - GSmt^ 

y +y H \ \ = (5.20) 

^ ^ A{t^-iy A{mt^-iy 4(t2-l)(mt2-l) ^ W uj 

Note that for this case, a + b = 4 and a — 6 = 3 are even and odd respectively. 
Hence, for such cases, one needs to use sets 2 and 4 from table 5.2 and sets 
1 and 3 from table 5.3, in order to get the four groups of the eigenfunctions. 
The solutions for this potential are given in table 5.5. 

We see that there is a degeneracy in the band edge energy eigenvalue 
14 — 7m + ^9 and all the solutions agree with the known solutions jS]. In 
this study, we have demonstrated the applicability of QHJ formalism to QES 
periodic potentials. We have been successful in obtaining the quasi - exact 
solvability condition and band edge solutions for cases 1 and 2. Case 3, 
requires some additional work and is not discussed here. 



Table 5.5: The residues, the value of n, number of linearly independent solutions, the band edge eigenfunc- 
tions and eigenvalues are as follows, with 5 — -\/4 — 4m + 25m^. Here, a + b — 4: and a — b — S, which give 
AT = 2 and M' = 1. 



set 




di 


n 


LI solutions 


Eigenfunction 


Eigenvalues 


1 


3/4 


1 


1 


1 


cnx(dn x)^''^snx 


5g — 171 + 2 


2 


3/4 





3 


2 


cna;(dn x)"^/^sn,T 
cna;(dna;)^"'^/^sna;(l — 2sn^a;) 


^9 — m + 2 
14 - 7m + ^9 


3 


1/4 


1 


3 


2 


(dna;)^/^(12msn^x - 5m - 2 - ^g) 
(dna;)^/2(12msn2x - 5m - 2 + ^g) 




2(5g 


4 


1/4 





4 


3 


(dnx)^/^(12msn "^x — 5m — 2 — 59) 
l^dnxf/^{12msn'^x -5m -2 + 69) 
1 — 8sn ^xcn '^x 




2(5g 

14 - 7m + (5g 
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Comparing the two different models of periodic potentials, namely the 
ES and QES, we find that the singularity structure of the QMF in both 
the cases is similar. For a potential with fixed values of a and b, these 
four conditions give four groups of solutions as seen as in table 5.2 and 5.3. 
A comparative analysis of periodic QES models with their non - periodic 
counterparts, reveals a few differences. For ordinary QES models, one obtains 
only one QES condition P in terms of n, which appears in the quantization 
condition ()2.70|) . Specifying the value of n, allows one to pick a potential 
from a family of potentials. Whereas, in the case of the periodic potentials, 
in general, one obtains four conditions as given in table 5.1. These four 
conditions can also be viewed as a constraint between a and b. If we select 
the value of n and the parameter a, one obtains four different values of b. 
For example, we choose the value a = 7/2 and n = A and substituting 
these values in the 4 conditions in table 5.1, one obtains the values of b 
as —5/2, 3/2, —3/2, 1/2. Among these, the values —5/2, 3/2 correspond to 
q = 15/4 and —3/2, 1/2, correspond q = 3/4. Thus, we have two potentials 
corresponding to distinct (p, g) values, namely (63/4,15/4) and (63/4,3/4). 
The form of the solutions for these potentials can be obtained from tables 
5.2 and 5.3. Bot these potentials will have a group of levels which have a 
QMF with four moving poles. 
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Chapter 6 



PT symmetric Hamiltonians 
6.1 Introduction 

Complex Hamiltonians possessing real eigenvalues have attracted consider- 
able attention in the current literature P-IH]. These quantal systems are not 
well understood because of their recent origin and apart from being counter 
intuitive. These Hamiltonians are characterized by a PT symmetry, which 
is discrete parity (x —x) followed by time reversal {i — > —i) symmetry. 
In the case, when the wave functions are also PT symmetric, the eigenvalues 
are real and the violation of PT symmetry by the wave function leads to 
eigenvalues which are complex conjugate pairs. Apart from identifying new 
Hamiltonians belonging to this class, the role of various discrete symmetries 
is also under thorough investigation. 

The presence of complex potentials in these systems, makes them ideal 
candidates to be probed using the QHJ formalism, since this approach has 
been formulated in the complex domain [TUl [TT| . It is extremely interesting 
to investigate the properties of the QMF of the PT symmetric Hamiltonians 

95 
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and to find out their diff'erences and similarities with ES real potentials and 
also the QES ones. 

In this chapter, we investigate the structure of the QMF of a class of 
PT symmetric Hamiltonians consisting of ES and QES models, for which 
the eigenfunctions and the eigenvalues are simultaneously obtained. The 
differences and similarities of these novel systems, with their ordinary ES 
and QES counter parts are clearly brought out. 

Application of the QHJ formalism to the PT symmetric potentials re- 
quires an approach different from that used earlier. In the absence of a 
generalization of the oscillation theorem, it is not clear whether the quanti- 
zation rule ()2.40|) is valid. If so what contour should be used. In this case, the 
quantization condition still holds for a contour which encloses all the moving 
poles. We shall assume this to be the case for the PT symmetric models to 
be taken up in this chapter. Here we will analyze two PT symmetric poten- 
tials namely, the Khare - Mandal potential in §6.2 and the complex Scarf 
potential V{x) = —Asech'^x — iBtanhx sechx, where A > in §6.3. 

6.2 Khare - Mandal model 

The potential expression for the Khare- Mandal model is V{x) = — (^cosh 2a;— 
iM^. This potential has either complex or real eigenvalues depending on 
whether M is odd or even |21 E] . Using the QHJ formalism, we obtain the 
QES condition, when M is odd and M is even. We also obtain the explicit 
expressions for the eigenvalues and eigenfunctions for the cases M = 3 and 
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M = 2. The QHJ equation in terms of g = ip, keeping h = 2m = 1 and after 
the change of variable t = cosh 2x becomes 

-dq 



+ 2V^^^ + E + {(t - tM)^ = 0. (6.1) 
The transformation equations in ()2.70|) for this change of variable become 

g = 2(v/t^)0, ^ = X- ^^J_^y (6.2) 
which transforms ()6.H) into 

^+rf^ + 4(^^+ 4(t2-l) =°' 
which will be treated as the the QHJ equation and x the QMF.. 

6.2.1 Form of QMF (x) 

We shall assume that x has a finite number of moving poles in the complex 
t plane and that the point at infinity is an isolated singularity. Besides the 
moving poles, x has fixed poles at t = ±1. It is seen from ()6.3j) that the 
function x is bounded at t = oo. Assuming that x has only these above 
mentioned singularities, we separate the singular part of x ci^nd write it in 
the following form. 

where hi and h'^ are the residues at fixed poles t = ±1 and -P„(t) is a polyno- 
mial of degree n. C gives the analytic part of x ^^id is a constant due to the 
Louville's theorem. From ()6.3|) . one can sec that x goes as ±i^/2, for large 
t, which are the values of C. 
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To find the residues at the fixed poles t = ±1, we expand x i^i Laurent 
series around each pole and proceeding in the same way as in the earlier 
chapters, one gets the following two values for bi and b[ as 

^ 3 1 3 1 

It has been assumed that the point at infinity is an isolated singularity. In 
order to find the residue of x infinity, one expands x terms of the 
Laurent series for large t as follows, 

X = ao + y + ^ + (6.6) 

and substitution of which in ()6.3p . gives the residue at oo as 

imC 

Ai = 6.7 
4ao 

along with oq = ±y, due to which the leading behaviour of x at infinity 
takes two values, which turns out to be 

A, = — . (6.8) 

Hence, using the fact, that for a rational function, the sum of all the residues 
is equal to zero, one obtains 

6i + 6' + n = Ai (6.9) 

From we see that the right hand side of ()(i.9j) is positive. Hence, for 

()6.9p to be true, we choose only the positive value of Ai = which fixes the 
value of residue at infinity. This means, we choose C = = +^^/2, since ao 
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gives the leading behaviour of x at infinity. It should be noted that there is 
no way of choosing a particular value of residue at the fixed poles, since one 
does not have information regarding the square integrability of the solutions. 
Hence, one needs to consider both the values of bi and b'l- Thus taking all 
possible combinations of bi and b[ in ()6.9|) . one obtains the QES condition for 
each combination along with a constraint on M as given in table 6.1. From 

Table 6.1: The QES condition and the number of moving poles of x for each 
combination of bi and b'l. 



set 


bi 


di 


n = X - 


- 6i - b\ 


Condition on M 


QES condition 


1 


1 

4 


1 

4 


M 
2 


1 

2 


M = odd, M >1 


M = 2n + l 


2 


3 
4 


3 
4 


M 
2 


3 
2 


M = odd, M > 3 


M = 2n + 3 


3 


3 
4 


1 
4 


M 
2 


- 1 


M = even, M >2 


M = 2n + 2 


4 


1 
4 


3 
4 


M 
2 


- 1 


M = even, M >2 


M = 2n + 2 



the table 6.1, we see that sets 1 and 2 are valid, only when M is odd and 
sets 3 and 4 are valid, only when M is even. 

6.2.2 Forms of the wave function 

Doing the change of variable and writing p in terms of x in ()2.74|) . one gets 
the equation for ipi^t) as 
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Hence, one can substitute the values of bi and b[ from sets 1 and 2 in (jfj.lOj) . if 
M is odd and sets 3 and 4 if M is even, to obtain the form of the wavefunction. 
The expression for the wave function is in terms of the unknown polynomial 
Pn{t), where n gives the number of zeros of Pn{t). In order to obtain the 
polynomial, we substitute x from ()6.4|) in ()6.3p to get 



2 



Pn{t) Pn{t) \t-l t+1 2 J {t-iy {t+iy 4(t2-i; 

E+iCt- iMf + 86i6; - AC\t^ - 1) + l±_ + A\ = 0. (6.11) 



4(^2-1) ^\t-l t+l^ 

This turns out to be equivalent to n linear homogeneous equations, for the 
coefficients of different powers of t in Pn{t). The energy eigenvalues are 
obtained by setting the corresponding determinant equal to zero. The explicit 
eigenvalues and eigenfunctions are obtained for M = 3 and M = 2 cases. 

6.2.3 Case 1 : M = 3 

Here, M is odd, so we can use sets 1 and 2 from table 6.1 and get the required 
results as below. 

Set 1 : bi = \, b[ = ^ and n = 1. 

This imphes P„ is a first degree polynomial, say, Bt + C. Substituting these 
values in 1)6.111) and comparing various powers of t, one obtains a 2 x 2 matrix 
for B and C as follows 



1 + ^ -< \(B 



0. (6.12) 
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Equating the determinant of this matrix to zero, one obtains the two values 
for energy and the polynomials as 

E = 7~e±2^1-Ae, Pi = |(2t-^(1±^1-4C2)). (6.13) 

Substituting the values of bi, b[ and Pi in (jfj.lOj) . gives the two eigenfunctions 
corresponding to the two eigenvalues as 

^fj{x) = e¥=°^^2a; |^2cosh2a; - ^(1 ± ^1 - 4(2) j . (6.14) 
Set 2: bi = |, b[ = | and n = 0. 

Here, we see that n = implies Pn(t) is a constant. Substituting these values 
in ()6.1H1 and proceeding in the same as before, one obtains 

E = 5 - C^ i^{x) = e^^°''^2^sinh2x, (6.15) 

which match with the known results |21ll]- 
6.2.4 Case 2 : M=2 

In this case, one makes use of sets 3 and 4 in table 6.1 and proceed in the 

same way as was done for case 1. 

Set 3 : bi = |, b\ = | and n = 0. 

The eigenvalues and eigenfunctions obtained are 

E = 3-C^ + 2iC, ^{x) =e'^^°'^^''{cosh2x+lY/\ (6.16) 

Similarly, for 

Set 4 ■' bi = |, b[ = I and n = 0, 
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one obtains 

E = 3-C^-2<, = e^™'^2^(cosh2x- 1)1/2. ^g^^^^. 

These match with the solutions given in 0]. Thus, for any given positive 
value of M, odd or even, one can obtain the eigenvalues and eigenfunctions 
for the Khare - Mandal potential. In the next section, we go over to the 
study of complex Scarf -II potential. 

6.3 Complex Scarf - II potential 

The expression for the potential is 

l/_ = — Asech^a; — ii?sech xtanhx. (6.18) 

The QHJ equation in terms of q, where g = ^ is 

+ - — V E + Asech x + iSsech xtanh x = 0. (6. 19) 

dx 

Doing the change of variable y = isinhx and proceeding in the same as 
was done for other the potentials, one obtains the QHJ equation for x, with 



F{y) = — y'^, as follows 

2 dx 2 + y^ _ E _ A + By ^ 
^ dy 4{l-y^y 1-y^ {1 - y^y ' ^' ' 

where the transformation equations ()2.70p are 
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Along with the n moving poles with residue one, x has poles at y = ±1. 
We assume that except for these poles, there are no other singularities and 
the point at infinity is an isolated singularity. Proceeding in the same way 
as in previous section, we obtain the residues at the poles and at infinity as 
follows. The residue at y = —1 



The residue at y = 1 



&i = ^±W^ + ^ + 5 (6.23) 



and the residue at infinity is 



A = i ± iVe. 



2 (6-24) 

As seen in the earlier section, one can write x in terms of its analytic and 
singular parts as 

h h' P' 

X = ^ + 4t + ^ + ^ (6.25) 
y-1 y+l Pn 

where, C is an analytic part of x ^"^^ is a constant due to Louville's theorem 
and it turns out to be zero. Using the fact, that the sum of all the residues is 
zero for a rational function, we obtain the expression for the energy eigenvalue 
as 

-E={hi + h\ + n-]^f. (6.26) 
The wave function in terms of x can be written using ()2.74|) as 



Chapter 6. PT symmetric Hamiltonians 104 
which is equal to 

^Ijiy) = iy-lf^--^iy+lf'^--^PM- (6-28) 

To obtain the polynomial Pn{y), one needs to solve, ()6.20|) by substituting x 
from ()(i.25|l . which gives the following differential equation 

P: + f-^ + + G(y)P„ = 0, (6.29) 

\y-l y+lj 

where 

^ im -h + h'^ - b[)) + 1 + AE + 8hb[)y' 
2y{^{h\ -h- bf + b[) - 2B)) 

(^2 _ 1)2 

{Ajpl -bi + b'^ - b[) + 2-AA-AE- 8bib[) 
+ (^^Tp • 

Substituting the expression for E from ()6.26p in ()6.29|) gives, 

(l-y2)P:+(2(6;-6i)-2(6i + 6;)y)P^+n(n+2(6i + 6;-l) + l)Pn = (6.30) 

which is in the form of the Jacobi differential equation and hence the poly- 
nomial Pn{y) is proportional to the Jacobi polynomial Pn^^ ^'^''^ ^{v)- Thus, 
the complete expression for the wave function in terms of the residues and 
the Jacobi polynomial in the x variable becomes 

ipi^x) = (isinhx - l)''i-^(isinhx + lf^-^P^^^'^-^'^^^-\ismhx). (6.31) 

Note that in this whole process, we had written the expression of the eigenval- 
ues and the eigenfunctions in terms of bi and b[, which have two values. No 
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particular value has been chosen. Hence, we need to choose one value of each 
residue to remove this ambiguity. For this case, unlike the Khare-Mandal 
potential, the solutions are known to be square integrable. Therefore, we 
make use of square integrability to choose one of the values of residues. We 
consider two different ranges of the potential parameters A and B and for 
each range, obtain the solutions. 

6.3.1 Case 1 : \B\> A+\ 

With this restriction on A and B, the residues at ?/ = ±1 becomes b[ = 
I ± ^\Jb — a — i and 6i = I ± A + B + ^ and thus, the wave function 
in ()f).31|l becomes 

iplx) = (isinhx — l)2^^(isinha; + 1)^^^ P^'^^'^'^^{ismh.x), (6.32) 

where r = B — A — j and s = ^ A + B + ^. From the above equation, one 
can see that for x — > cxd, ?/^(a;) ^ 0, for 

with n limited to lie between < n < A + B + ^ — ^. Thus, for this 
parameter range, expression for tplx) will be 

ip{x) = {isinhx — l)3"^(isinha; + l)3=^^P^*'''~''(isinha;) (6.34) 

and the corresponding expression for energy from ()6.26|) will be 

E = -(n+l~\(^|^7^- + ±^^^^I^]] . (6.35) 
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6.3.2 Case 2 : \B\<A + \ 

In this case, the wave function in (jfj.Hlll is square integrable, if 

''^-2-2^^^^-"' 2- 2^^ '^'■''^ 

The expression for the wave function is given by, 

iplx) = (isinhx — 1)4^2 (^isinhx + l)4~2p~* "''(^sinha:), (6.37) 

where t = ^J^ + A — B and s = ^J^ + A + B and the energy eigenvalue 
expression in ()6.26|1 becomes 

with n restricted by 

0<n< y^ + A-B + yA + B + ^ - i (6.39) 

In conclusion, PT symmetric potentials belonging to the QES and ES 
class have been investigated through the QHJ formalism. The QES solv- 
able Khare-Mandal potential has complex or real eigenvalues, depending on 
whether the potential parameter M is odd or even. The singularity structure 
of the QMF for these two cases is different. For the case, when M is odd, one 
observes from table 6.1, that the solutions fall into two groups, which consist 
of solutions coming from sets 1 and 2. For a solution belonging to a partic- 
ular group, the number of singularities of the QMF are fixed and consists of 
both real and complex locations. This kind of singularity structure of the 
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QMF has been observed in the study of periodic potentials [121 CBl- Though 
the solutions, for M even, fall into two groups coming from sets 3 and 4, 
they all have the same number of singularities, which again can consist of 
complex and real poles. This singularity structure is same as that observed 
in the ordinary QES models [T^ . 

Coming to the case of exactly solvable PT symmetric potential, the lo- 
cation of the moving poles can be either real or complex. In the specific 
example of complex Scarf potential, it turns out that all the moving poles 
are off the real line. In contrast for the ordinary ES models fS] , the moving 
poles are always real. For both the cases, the number of moving poles of the 
QMF characterize the energy eigenvalues. 
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Chapter 7 



Conclusions 



In this thesis, we have studied various potential models with special spectral 
properties using QHJ formalism. The various models, being potentials which 
exhibit different spectra for different ranges of parameters, potentials with 
band structure and potentials with both real and complex eigenvalues. For 
all the models studied, we obtained the eigenvalues and the eigenfunctions 
successfully. Thus, the QHJ formahsm provides an alternate method to solve 
various types of ES and QES models. The added advantage of this method is, 
that it is very simple and makes use of well known complex variable theorems. 

The most important steps in obtaining the solutions, for both ES and 
QES models, have been, the choice of the change of variable and use of 
the singularity structure of the QMF in the complex plane. We did not 
attempt to derive the singularity structure of the QMF within the QHJ 
formalism. However, we assumed that the singularity structure is very simple 
viz. the QMF has a finite number of moving poles. With this assumption, 
we could obtain the bound state and band edge solutions in a most straight 
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forward fashion, even if the equations to be solved were not of a simple form. 
The fact that one arrived at correct solutions in all the cases, supports a 
conjecture that the assumption may be true for all ES and QES models. It 
should be remarked that this assumption is equivalent to the property that 
'the QMF', becomes a rational function after the suitable transformations. 
Stated differently, this assumption is equivalent to, the point at infinity being 
an isolated singular point. In this work our study has been limited to one 
dimensional models. It will be interesting to investigate if any of these forms 
of the assumption remain valid for higher dimensional models. 

In all the models, the integer n appearing in the exact quantization con- 
dition is the number of moving poles of the QMF. For ES models P and 
for Scarf potential, these poles are located on the real axis only, with no 
moving poles off the real axis. Given a value of n, there is a unique energy 
level and the corresponding eigenfunction has only real, n, zeros. However, 
for the QES models |2], n appears as a parameter in the potential and of 
the infinite possible states, only a finite number, determined by n, of states 
can be obtained analytically. In the case of QES models, the moving poles 
appear at complex locations also. Specifying a value of n, selects a potential 
within a family and does not pick up a unique level for the potential. In fact, 
all the analytically solvable eigenfunctions have the same number (= n) of 
zeros, of which the number of real zeros are in accordance with the oscillation 
theorem [H]. 

The periodic potentials in chapters IV 0], however, show a richer distri- 
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bution of the moving poles in the complex plane. Like ordinary QES models, 
the moving poles of the QMF are located both on and off the real axis. Here, 
in this case, selecting a value of the integer n, picks up a subset of solutions 
all having n zeros just like QES models. However, unlike the QES models, 
varying the integer n does not give a new potential, instead it gives a different 
subset of levels. 

The case of QES periodic potentials jH] of chapter V is very similar to 
the ES periodic potentials, as far as the singularity structure of the QMF is 
concerned. 

The above conclusions have been arrived in the context of specific models 
studied. However, we expect them to remain generally valid for other one 
dimensional, periodic and aperiodic, ES and QES potentials. The number 
of PT symmetric models [0] studied here is too small to allow any similar 
generalizations. 

Our study exhausts analytically solvable bound state and band edge solu- 
tions for both ES and QES cases of one dimensional potentials. For problems 
which are not ES, making a proper guess about the singularity structure of 
the QMF, does not appear to be possible. A study of the distribution of 
the moving poles in the complex plane may have a relation to the classical 
properties of the system. Such a study is interesting in its own right and 
can be best done numerically. We expect that, here also, the established 
results in the complex variable theory will provide a useful scheme to obtain 
numerical, possibly approximate analytical solutions for the wave functions 
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and energies. Finally, we hope that a useful extension to continuous energy 
states in one dimension and to non separable models in higher dimensions 
can be found. 
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